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O ' Abstract 

Using methods of effective field theory, factorized expressions for arbitrary B — > X u I v 
decay distributions in the shape-function region of large hadronic energy and moderate 
hadronic invariant mass are derived. Large logarithms are resummed at next-to-leading 
order in renormalization-group improved perturbation theory. The operator product 
expansion is employed to relate moments of the renormalized shape function with HQET 
parameters such as m&, A and Ai defined in a new physical subtraction scheme. An 
analytic expression for the asymptotic behavior of the shape function is obtained, which 
reveals that it is not positive definite. Explicit expressions are presented for the charged- 
lepton energy spectrum, the hadronic invariant mass distribution, and the spectrum in 
the hadronic light-cone momentum P + = En — \Ph\- A new method for a precision 
measurement of \V u b\ is proposed, which combines good theoretical control with high 
efficiency and a powerful discrimination against charm background. 



1 Introduction 



Inclusive _B-meson decays into light final-state particles, such as the semileptonic process 
B — > X u l~v and the radiative process B — > X^, are of great importance to the determination 
of the Cabbibo-Kobayashi-Maskawa (CKM) matrix element |T4&| and play a prominent role in 
the search for New Physics at the B factories. The ever increasing accuracy of experimental 
data on the inclusive decay rates and spectra provides a strong incentive to continuously 
improving the theoretical description of these processes. To keep up with the reduction of 
experimental errors, it will soon be necessary to have theoretical predictions for inclusive 
decay distributions with uncertainties below the 10% level. 

There are several challenges facing theorists pursuing this ambitious goal. Due to ex- 
perimental cuts one is generally faced with a situation where the hadronic final states arc 
constrained to have large energy E H ~ M B but only moderate invariant mass s# ~ M b Aq CD . 
In this kinematic region (called the "shape-function region"), the decay rates and spectra are 
obtained using a twist expansion [1, 2, 3], which is considerably more complicated than the con- 
ventional heavy-quark expansion employed in the calculation of inclusive rates for B — > X c l~v 
decays [4, 5, 6]. In the twist expansion, infinite sets of power corrections are resummed into 
non-perturbative "shape functions" describing the internal structure of the B meson. At the 
same time, the short-distance corrections are also more complicated than in most other ap- 
plications of heavy-quark expansions, because the kinematics implies the relevance of three 
separated mass scales: M B ("hard"), -J Mb Aqcd ("hard-collinear"), and Aqcd ("soft"). Large 
logarithms of ratios of these scales arise at every order in perturbation theory and must be 
resummed [7, 8, 9]. To properly disentangle the physics associated with these scales requires 
a sophisticated effective field-theory machinery [10], which has only been fully developed re- 
cently. A systematic treatment consists of matching QCD onto soft-collinear effective theory 
(SCET) in a first step, in which hard quantum fluctuations are integrated out. In a second 
step, SCET is matched onto heavy-quark effective theory (HQET), and hard-collinear modes 
are integrated out. The resulting expressions for inclusive differential decay rates have the 
factorized form dY ~ H J <g) S [8, 9]. The function H contains the hard corrections, the jet 
function J, which describes the properties of the final-state hadronic jet, contains the hard- 
collinear effects, and the shape function S accounts for the internal soft dynamics in the B 
meson [1, 2]. The <g> symbol implies a convolution over a light-cone momentum variable u 
associated with the residual momentum of the b quark inside the B meson. 

Up until now, no complete next-to-leading order predictions for inclusive decay rates and 
spectra in the shape-function region have been presented in the literature. In the present 
paper we close this gap. We go beyond previous work in several important ways. First, 
we complete the matching calculations for the two-step matching QCD — > SCET — > HQET at 
next-to-leading order in perturbation theory. We then derive renormalizat ion-group equations 
governing the dependence of the functions H, J, and S on the renormalization scale, and solve 
these equations analytically in momentum space. Next, we derive several model-independent 
properties of the shape function S, which were so far unknown. In particular, we present 
the precise form of the relations between renormalized shape-function moments and HQET 
parameters such as A and Ai, and we give an analytical formula for the asymptotic behavior 
of the shape function. An unexpected outcome of this analysis is the finding that the shape 
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function is not positive definite, but acquires a negative radiative tail at large values of 
We also discuss how to derive shape-function independent relations between different decay 
spectra, which are free of spurious Landau-pole singularities. This improves on similar relations 
that can be found in the literature [11, 12, 13]. Finally, we propose a new method for a high- 
precision determination of \V u b\, which offers several advantages over previous approaches for 
extracting this important CKM parameter. 

The remainder of this paper is organized as follows: In Section 2, we present a derivation 
of the factorization formula for inclusive rates in the shape-function region using the position- 
space formulation of SCET. The one-loop matching calculations needed to derive perturbative 
expressions for the functions H, J, and S are presented in Section 3. In Section 4, we derive 
evolution equations for these functions and present their exact analytic solutions in momen- 
tum space. A comprehensive discussion of the properties of the shape function is given in 
Section 5. While much of the discussion in this paper is necessarily rather technical, our 
results have important implications for the phenomenology of inclusive 5-meson decays, in- 
cluding in particular the determination of \V u b\. We have therefore organized the paper in 
such a way that the reader interested mainly in the applications of our formalism can skip the 
conceptual parts of the discussion in Sections 2-5 and proceed directly with Section 6, where 
we summarize our findings and collect all formulae needed for phenomenology. The following 
Sections 7, 8, and 9 contain explicit predictions for various decay rates and spectra, which are 
of relevance to experimenters. 

During the final stages of this work, a paper by Bauer and Manohar appeared [14], which 
overlaps with parts of our analysis. These authors derive expressions for the hard-scattering 
kernels, the jet function, and the renormalized shape function which agree with our results 
reported in (25), (27), and (33) below. They also obtain an expression for the anomalous 
dimension of the shape function, which coincides with our result in (43). In many other 
aspects our analysis goes beyond that of [14]. In particular, using a technique published some 
time ago by two of us [15], we succeed to obtain an exact analytic solution of the evolution 
equation for the shape function in momentum space. At leading-logarithmic order we confirm 
an expression found earlier by Balzereit, Mannel, and Kilian (apart from a small mistake) [16]. 
While in [14] resummed expressions are only given for large- N moments of decay spectra, we 
present results for arbitrary decay distributions in physical phase space. We also disagree 
with Bauer and Manohar on the conclusion that moments of the renormalized shape function 
cannot be related to HQET parameters. In fact, we derive the explicit form of such relations. 

2 Factorization theorem for inclusive decay rates 

Using the optical theorem, the hadronic physics relevant to the inclusive semileptonic decay 
B — > X u l~v can be related to a hadronic tensor W^ v defined via the discontinuity of the 
forward 5-meson matrix element of a correlator of two flavor- changing weak currents J M = 
uY(l - 75)6 [4, 5, 6]. We define 
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Here v is the 5-meson velocity and q the momentum carried by the lepton pair. The current 
correlator receives contributions from widely separated energy and distance scales. To extract 
the dependence on the large 6-quark mass, it is convenient to rescale the heavy-quark field 
according to b(x) = e~ imbV ' x b' (x) . The field b'(x) carries the residual momentum k = pb — 
m b v = O(Aqcd)- Then the phase factor in (1) becomes e l (<i-m b vyx = e -ip-x j n ^ e p ar ^ on 
model, p = rribV — q corresponds to the momentum of the jet of light partons into which the 
6-quark decays. 

As long as at least some components of the jet momentum are large compared with Aqcd, 
the current correlator can be evaluated using a short-distance expansion. The simplest case 
arises if all components of p^ ~ m h ^> Aq CD . Integrating out the light-quark propagator then 
produces bilocal operators of the form b'(x) . . . fe'(O) with x ~ 1/mb, which can be expanded in 
terms of local operators O n multiplied by coefficient functions D n (mb, v-x, x 2 ). Performing the 
integration over x gives functions D n (rrib, v ■ p, p 2 ) = J d 4 x e~ ip ' x D n (rrib, v ■ x, x 2 ) depending on 
the different short-distance scales [5, 6]. The 5-meson matrix elements of the local operators 
O n can be evaluated using techniques of HQET [17]. Thus, the operator product expansion 
produces an expansion in inverse powers of the heavy-quark mass. This expansion is relevant 
to the computation of inclusive decay rates without restrictive cuts on kinematic variables. 
Such rates are governed by parton kinematics as imprinted through the dependence on the jet 
momentum p. The unique operator of lowest dimension is hh (with h the effective heavy-quark 
field in HQET), whose forward matrix element between 5-meson states is unity. The Wilson 
coefficient of this operator has been computed to 0(a s ) in [18]. At next-to-leading order in the 
heavy-quark expansion two new operators arise, which correspond to the kinetic energy and 
chromo- magnetic interaction of a heavy quark. Their Wilson coefficients have been computed 
at tree level in [5, 6]. 

A more complicated situation arises in the region of phase space in which the energy of the 
hadronic final states is much larger than its invariant mass, meaning that some components of 
the vector p^ are much larger than others. In this case, the current correlator can be expanded 
in non-local light-cone operators [1, 2, 3]. The formalism is most transparent when presented 
using the language of soft-collinear effective theory (SCET) [10]. To simplify the kinematics 
it is convenient to work in the 5-meson rest frame, where t> M = (1,0,0,0), and to choose 
the lepton momentum q along the negative z-axis, so that p points in the z direction. Next, 
we define two light-like vectors n M = (1,0,0,1) and n M = (1,0,0,-1) satisfying n • n — 2, 
n • v — n • v — 1. Any 4- vector can be expanded in the light-cone basis as 

= (n ■ p) — + (n ■ p) — + = p\ + j£ + p%_, (2) 

where p± ■ n = p± ■ n = 0. In the shape-function region, the jet momentum (and likewise the 
momentum of the hadronic final state) scales like p^ ~ E(X,1,\/X), where A ~ Aqqd/E is 
the SCET expansion parameter and E ~ rrib ^> Aqcd- (For the jet momentum, p± = by 
choice of the coordinate system.) It follows that pt ~ E is a hard scale whereas p\ ~ Aq CD 
is a long-distance hadronic scale. The jet invariant mass, p 2 ~ -EAq CD , defines a hybrid, 
intermediate short-distance scale. We refer to a momentum with these scaling properties as 
"hard-collinear" . The appropriate effective field theory for integrating out the short- distance 
fluctuations associated with the hard scale p_ is called SCETi and has been discussed in 
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detail in the literature [19, 20, 21, 22]. Here we use the coordinate-space formulation of SCET 
developed in [22, 23]. 

Below a matching scale /^ ~ m b , the semileptonic current can be expanded as 

3 r ~ 

u(x)^(l -is)b'{x) = dsC l (s)X(x + sn)T t ;K(x„) + ... , (3) 

i=i J 

where the dots denote higher-order terms in the SCET expansion, which can be neglected 
at leading power in Aq C d/^6- The soft heavy-quark field Ji(x-) = Sl(x-) h(x ) and the 
hard-collinear light-quark field X(x) = Sl(x_) Wl c {x)t;{x) are SCET building blocks that are 
invariant under a set of homogeneous soft and hard-collinear gauge transformations [23, 24, 25]. 
The objects S s and Wh c are soft and hard-collinear Wilson lines [20]. Soft fields in SCET 
are multi-pole expanded and only depend on the minus component of the position vector 
x. Finally, the position-space Wilson coefficient functions Cj(s) depend on the variable s 
defining the position of the hard-collinear field. Here and below we denote functions in position 
space with a tilde, which is omitted from the corresponding Fourier-transformed functions in 
momentum space. A convenient basis of Dirac structures in (3) is 

1^ = 7^(1-75), 1^ = ^(1 + 75), It = — (l + 7fe)- (4) 
The current correlator in (1) then becomes 

T^ = i J d 4 xe~ ip - x J2 J dsdtC*(t)C i (s)T{M(0)r»X(tn),X(x + sn)r" i K(x-)} + ... . 

(5) 

In a second step, the hard-collinear fluctuations associated with the light-quark jet can be 
integrated out by matching SCET onto HQET at an intermediate scale fa ~ y'm^AQCD. At 
leading order the SCET Lagrangian (when written in terms of the gauge-invariant fields such 
as X) does not contain interactions between hard-collinear and soft fields. Since the external 
S-meson states only contain soft constituents, we can take the vacuum matrix element over 
the hard-collinear fields, defining a jet function 

(Q\T {X k (tn),Xi{x + sn)} \Q) = 8 kl J(x + (s - t)n) + . . . , (6) 

where k, I are color indices, and we have used translational invariance to determine the de- 
pendence on the coordinate vectors. At higher orders in SCET power counting, additional 
jet functions would arise, but their contributions can be neglected at leading power. Shift- 
ing the integration variable from x to z = x + (s — t)n, with Z- — X-, and introducing the 
Fourier-transformed Wilson coefficient functions 

Ci(n-p) = j dse ian - p Ci(s), (7) 

we then obtain 

3 r 
T^ = t J2 C* (n ■ p) Q (n-p) / d A z e"^ ft(0) f£ J{z)T\Ji{z-) + .... (8) 

i,3=l J 
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In the next step, we rewrite the bilocal heavy-quark operator as [1] 

X(0)rK(z_) = (h S S ){0)T e z - d + {St h){0) = h{0)T e z - D + h{0) 

diue-^ n - z h(0)T5(cu-in- D)h(0), (9) 



where T may be an arbitrary (even dependent) Dirac structure, and we have used the 
property in ■ D S s = S s in • d of the soft Wilson line S s , where iD^ = id^ + g s A^ is the 
covariant derivative with respect to soft gauge transformations. When this expression is used 
in (8), the resulting formula for the correlator involves the Fourier transform of the jet function, 



/ 



d 4 ze- ip - z J(z)=^J(p 2 ), (10) 



however with p^ replaced by the combination p^ = p^ + ^cun^. In defining the momentum- 



space jet function J{p 2 ) we have taken into account that the matrix element in (6) vanishes 
when multiplied by ifi from either side, since yiX = 0. Using that — p_, we now obtain 

3 

T"" = iJ2 H a(n-P) / duJ{pl)hT^_T v i 8{u-in-D)h + ... , (11) 

where Hij(n ■ p) = C*{n • p) Ci{n • p) are called the hard functions. 

In order to compute the hadronic tensor we take the discontinuity of the jet function, 

J(p 2 ) = - Im [tj(p 2 )] , (12) 

7T 

and evaluate the 5-meson matrix element of the soft operator using the HQET trace formalism, 
which allows us to write [17] 

(g WI hT6(.-in.D)h |g M ) = 1 t 1+A 
2M B 2 V 2 / 

at leading power in the heavy-quark expansion. The soft function S(u) coincides with the 
shape function f(k+) introduced in [1]. This gives the factorization formula 

W" = H ^ ■ P) tr (f J t rr i±^) | J(pl) S(u>) + .... (14) 

At leading power one could replace $L — > f 1 and n ■ p = 2v ■ p_ — > 2t> • p in this result. 

In the final expressions (11) and (14) the dependence on the three scales n ■ p ~ m&, 
~ m 6 Aq CD and a; ~ A QCD has been factorized into the hard, jet, and shape functions, re- 
spectively. Large logarithms associated with ratios of these scales can be resummed by solving 
renormalization-group equations for the scale dependence of these component functions. The 
factorization formula (14) was derived at tree level in [1, 2], and was generalized to all orders 
in perturbation theory in [8, 9]. The derivation presented above is equivalent to a proof of this 
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formula presented in [20] (see also [10]). The limits of integration in the convolution integral 
are determined by the facts that the jet function defined in (12) has support for p^ > 0, and 
the shape function defined in (13) has support for — oo < uo < A with A = Mb — nit,, where 
rrib is the heavy-quark pole mass. The argument of the jet function can be rewritten as 

p 2 w = p 2 + n ■ puj = n ■ p(n ■ Ph — (A — uj)) = n ■ p(n ■ Ph — Cj) , (15) 

where Ph = Msv — q = p + Av is the 4- momentum of the hadronic final state, and the variable 
Cj = A — uj > 0. Finally, n-Pn = Eh — \Ph\ = s #/2i?ff + 0(AQ CD /m{,) is a kinematic variable of 
order Aqcd related to the hadronic invariant mass and energy of the final state. The usefulness 
of this variable has also been emphasized in [16, 26, 27]. We shall see below that expressing 
the convolution integral in terms of the new variable Cj eliminates any spurious dependence of 
the decay spectra on the 6-quark pole mass. It follows that in the shape-function region the 
hadronic tensor is most naturally written as a function of the "parton variable" n ■ p and the 
"hadron variable" n ■ Ph, and the inclusive spectra are governed by a combination of parton 
and hadronic kinematics. 

Using the fact that the Wilson coefficients Cj are real and hence is symmetric in its 
indices, we find 

Hij tr (i* t r- l±f \ = 2H n (p»_v» + p^ - v ■ p_ - ie^p^vp) (16) 
i,j=i ^ • ' 

+ 2H 22 v ■ p. v»v v + 2(ffi 2 + H 23 ) {f_v v + p v _v») + 2(2H 13 + H 33 ) . 

v ■ p- 

This result may be compared with the general Lorentz decomposition of the hadronic tensor 
given in [18]: 

W» v = W 1 {fv v + p v v» -g^vp- ie» vaP p a vp) - W 2 g^ 

+ W 3 v"v u + W 4 (p^v u + p u v") + W 5 p»p v (17) 

We see that the structure function W 2 is not generated at leading order in the SCET expan- 
sion. Since only the Wilson coefficient C\ is non-zero at tree-level, the structure function W\ 
receives leading-power contributions at tree level, whereas W± and W$ receive leading-power 
contributions at 0(a s (mi>)). The function W 3 receives leading-power contributions only at 
0(ct!g(ra&)), which is beyond the accuracy of a next-to-leading order calculation. 

3 Matching calculations 

In this section we derive perturbative expressions for the hard functions Hij(n ■ p) and the 
jet function J(p^) in (14) at next-to-leading order in a s . To this end, we match expressions 
for the hadronic tensor obtained in full QCD, SCET, and HQET, using for simplicity on-shell 
external 6-quark states. We also present results for the renormalized shape function in the 
parton model, which are needed in the matching calculation. Throughout this paper we use 
the MS subtraction scheme and work in d = 4 — 2e space-time dimensions. 
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3.1 Hard functions 



Perturbative expressions at 0(a s ) for the structure functions Wi in the decomposition (17) 
have been obtained in [18] by evaluating one-loop Feynman graphs for the current correlator 
T^ v using on-shell external quark states with residual momentum k (satisfying v ■ k — 0) in 
full QCD. The leading terms in the region of hard-collinear jet momenta are 



- Wi 
2 



Cpot 



C F a s 



in 



8 In 2 y - 10 my + + 4L 2 (1 - y) + ^ + 5 



+ 



47T 



ln(p 2 /r< 



^ 4 



47T 1 — y V 1 — y 



2 Cj?a s 



+ (81ny-7) 
f 1 ) + . . • , 



3 

+ ... , 



(18) 



m b 



47r 1 — y \ 1 — y 



l-2y 



lny-1 +... , 



whereas W2 and W3 do not receive leading-power contributions at this order, in accordance 
with our general observations made above. Here a s = a s (fi), y = n-p/mb, andpf, = p 2 +n-pn-k. 
We have used that 2v ■ p/rrib = n ■ p/rrib + 0(A) in the hard-collinear region. The star 
distributions are generalized plus distributions defined as [18] 



z , F(x)-F(0) „ /n .,z 
dx — + F(0) In - , 



x 



u 



dxF(x) 



<o 



ln(x/u)\^ u 



x 



'■ F(x)-F(0) , x F(0) , 2 2 
da; — In — I — In 2 - 

x u 2 u 



(19) 



where F(x) is a smooth test function. For later purposes, we note the useful identities 



A 



>(i) 



\n(Xx/u) 
\x 



[u/X] 



+ S(x) In A, 
\n(Xx/u) Y"/ A ' f\n(x/u 



x 



x 



+ 



In A + In 2 A 



(20) 



In order to find the hard functions if^, we calculate the discontinuity of the current 
correlator (5) between on-shell heavy-quark states in momentum space. We work to one- 
loop order in SCET, keeping i,j fixed and omitting the Wilson coefficient functions. The 
corresponding tree diagram yields 



£>(°) =K8(pl), with K = u b (v) f 1 f_ Y\ u b (v) , 



(21) 



where u b {v) are on-shell HQET spinors normalized to unity, and the quantity K corresponds 
to the Dirac trace in (14). The interpretation of this result in terms of hard, jet, and soft 
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Figure 1: One-loop diagrams contributing to the current correlator in SCET. 
The effective current operators are denoted by crossed circles, and hard- 
collinear propagators are drawn as dashed lines. Mirror graphs obtained by 
exchanging the two currents are not shown. 



functions is that, at tree level, J^\pt) = &{pt) and S^ rton {uo) = 5(uj—n-k). (The second result 
is specific to the free-quark decay picture.) Then the convolution integral J dou J{plj) S(uj) in 
(14) produces S{pl), and comparison with (18) shows that = 1, while all other hard 
functions vanish at tree level. 

The diagrams contributing at one-loop order are shown in Figure 1. They are evaluated us- 
ing the Feynman rules of SCET. The first three graphs contain hard-collinear gluon exchanges, 
while the last two diagrams contain soft exchanges. The wave-function renormalization fac- 
tors of the external heavy quarks equal 1 on-shell. For the sum of all hard-collinear exchange 
graphs, we find 



+ - + 7-7T 2 )6(pt)+4 



he 4n 



The sum of the soft contributions is given by 



pI 



[M 2 



+ 3 



pIJ* 



(22) 



= K 



C F a s 
An 



4 



2 

e 



4L 2 + AL 



MpI/v 2 

Pi 



In 2 



+ - + 8L - 4 



71 

~6 



Pi 



5(pl) 



[m 2 



(23) 



where L = ln(n-p//i). The 1/e poles in the sum of the hard-collinear and soft contributions are 
subtracted by a multiplicative renormalization factor Zj applied to the bare current correlator 
in (5), where 

C F a s ( 1 2 r 5 \ 



47T 



is the (momentum-space) current renormalization constant in SCET [19]. Taking the sum of 
(22) and (23) after subtraction of the pole terms, and matching it with the results in (18), we 
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find that at one-loop order 

H 11 (n-p) = l + ^p- (-AL 2 + 10L - 4 In y - ^ - 4L 2 (1 - y) - ^ - 12 
47r \ 1 — y o 

w i- \ Cf(Xs 2 ( y ln y 1 

^• p) = ^rr^iT^ + 1 J' (25) 

„ /_ \ C F a s y /l - 2y 

Hiz{n-p) = — lny- 1 

An 1-y \l-y 

In deriving these results we have used the identities (20) to rearrange the various star distri- 
butions. The remaining hard functions start at 0(a 2 ). Using the relation Hij = CiCj, one 
can derive from these results expressions for the Wilson coefficients in the expansion of the 
semileptonic current in (3). We confirm the expressions for these coefficients given in [19]. 



3.2 Jet function 



After the hadronic tensor is matched onto HQET as shown in (14), the SCET loop graphs 
in Figure 1 determine the one-loop contributions to the product of the jet function and the 
shape function in (14). We may write this product symbolically as <g> + ® 
where the <g> symbol means a convolution in uj. Whereas the jet function is a short- distance 
object that can be calculated in perturbation theory the shape function is defined in terms of 
a hadronic matrix element and cannot be properly described using Feynman diagrams with 
on-shell external quark states. In a second step, we must therefore extract from the results 
(22) and (23) the one-loop contribution to the jet function. To this end, we must compute 
the renormalized shape function at one-loop order in the parton model. This will be done in 
the next subsection. We may, however, already anticipate the result for the jet function at 
this point, because at one-loop order the graphs in Figure 1 can be separated into diagrams 
with hard-collinear (first line) or soft (second line) gluon exchange. (This separation would be 
non-trivial beyond one-loop order.) The hard-collinear contribution in (22) thus determines 
the convolution ® = whereas the soft contribution in (23) corresponds to 

® S^ 1 ). It follows that the renormalized jet function is given by the distribution 



J(PL) = 6(Pt) + 



CfOC s 

An 



(7 - vr 2 ) 5(pi) + 4 



vl 



-3 



(26) 



This result disagrees with a corresponding expression obtained in [28]. It will often be useful 
to separate the dependence on n ■ p and n ■ Ph in this result by means of the substitution 
vl — VVw-, where p^ — m^n • Ph — &) according to (15). Using the identities (20), we find 



yj(pl) = j(pl,y) = 6(pl) + 



C F a s 
An 



(2\n 2 y-3\ny + 7-n 2 )5(p 2 J 



+ 4 



fMvl/v 2 



i^ 2 



+ (41ny-3) 



y 2 



V vi A ' Wu 

The jet function is non-zero only if y > and n • Ph > &, which ensures that p 2 ^ > 0. 



(27) 
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Figure 2: Radiative corrections to the shape function. The bilocal HQET 
operator is denoted by the black square. A mirror copy of the first graph is 
not shown. 



3.3 Renormalized shape function 



Having calculated the short- distance objects Hij and J in the factorization formula (14) at 
one-loop order, we now turn to a study of radiative corrections to the shape function S(u). 
There is considerable confusion in the literature about the renormalization properties of the 
shape function, and several incorrect results for its anomalous dimension have been published. 
We therefore present our calculation in some detail in this subsection and the following section. 

According to (13), the shape function is defined in terms of a hadronic matrix element in 
HQET and thus cannot be calculated perturbatively. However, the renormalization properties 
of this function can be studied using perturbation theory. To this end, we evaluate the matrix 
element (13) in HQET using external heavy-quark states with residual momentum k. For 
the time being, we keep v ■ k non-zero to regularize infra-red singularities. The relevant one- 
loop graphs are depicted in Figure 2. Adding the tree contribution, we obtain for the matrix 
element of the bare shape-function operator 0{uj) = hT 5(u — in - D) h (expressed in terms of 
renormalized fields) 



Share M = Z h 5(u) - 71- k) 



4Cf9 2 s 
(4tt) 2 - 



r(i + e) 



x 



dll 



-l-2e 



5{ijj — n ■ k) — 5{ui — n ■ k + I) 



1 + 



+ 6{n • k — ui) (n ■ k — ui) e {n • k — uj + 5) 



-l-e 



(28) 



where 5 = —2v ■ k, and 



1 + A ° F f\ r(2e) T(l - e) 5~ 2e 



(29) 



(4tt) 2 - £ 

is the off-shell wave-function renormalization constant of a heavy quark in HQET. The next 
step is to extract the ultra-violet poles from this result, which determine the anomalous di- 
mension of the shape function. We define a renormalization factor through 



dJ Z s (u,uj') S ha , rc (uj') , 



C F a s 



(30) 



Z s (to, uj') = 5(cu - uj') H — - zs(uj, uj') + . . . 

4"7T 
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The result for Z s following from (28) must be interpreted as a distribution on test functions 
F(uj') with support on the interval — oo < uj' < A. We obtain 



z s (u,cf) 



-1)s(uj-uj')-- 
e 2 e A-uj e e 



2 4 , 
+ - In 



6{uj' - uj) 
uj' — UJ 



UJ 



(31) 



Note the peculiar dependence on the parameter A setting the upper limit on the integration 
over uj' in (30), which combines with the plus distribution to form a star distribution in the 
variable (u/ — uj). 

We can now determine the renormalized shape function from (30). The result must once 
again be interpreted as a distribution, this time on test functions F{uj) integrated over a finite 
interval — A had < uj < A. In practice, the value of A had is set by kinematics or by virtue of 
some experimental cut (see Sections 7 and 9 below). The result is 



Spartan M = 5(uJ - Tl ■ k) 



C F a s 



C F a s 



7Y 



7T 2 / -5 

h L 2 [ 

24 \A had + n-k 



7Y 



6{n ■ k — uj) 
n ■ k — uj 



In 



n ■ k 



uj , n • k 
h ln- 



UJ 



+ 5(n-k-uj) In 2 



Ahad + n ■ k 6{n • k — uj) 



// 



+ 



n ■ k — uj + 5 



+ 8(n ■ k — uj) In 



5 



(32) 



While it was useful to keep the heavy quark off-shell in the calculation of the ultra-violet 
renormalization factor, the limit 5 = —2v-k — > can be taken in the result for the renormalized 
shape functions without leading to infra-red singularities. This gives 



'parton 



5{uj — n • k) 
C F a s 



1 - 



C F a, 



7T 



7T 



7T 



n ■ k — uj 



In 



24 
n ■ k — uj 



+ 



n ■ k — uj 



(33) 



where the star distributions must now be understood as distributions in the variable (n-k — uj). 

We stress that these results for the renormalized shape function are obtained in the parton 
model and can in no way provide a realistic prediction for the functional form of S(uj). This 
should be obvious from the fact that our results depend on a single "hadronic parameter" n-k, 
corresponding to a fixed residual momentum of the heavy quark. Only the dependence on 
the ultra-violet renormalization scale /i can be trusted. However, the one-loop result in (33) 
is needed to complete the matching calculation of the jet function described in the previous 
subsection, which can legitimately be performed with on-shell external 6-quark states. Given 



11 



the expression for the renormalized shape function, we obtain 

/ dujJ (0 \pl) S (1 \lo) = —S {1) (-p 2 /n-p) (34) 
J n-p 



C F a s 



TV 



24 n-p V p{ n-pfij^ n-p\ p{ 



With the help of the identities (20) this can be shown to be equal to the finite part of (23), 
as we claimed above. 

4 Renormalization-group resummation 

Equations (25) and (26) determine the short-distance objects and J in the factorization 
formula (14) at one-loop order in perturbation theory. However, there is no common choice 
of the renormalization scale /i that would eliminate all large logarithms from these results. 
Likewise, the shape function, being a hadronic matrix element, is naturally renormalized at 
some low scale, whereas the short-distance objects contain physics at higher scales. The 
problem of large logarithms arising from the presence of disparate mass scales can be dealt 
with using renormalization-group equations. Proceeding in three steps, our strategy will be 
as follows: 

i) At a high scale Hh ~ n^b we match QCD onto SCET and extract matching conditions for 
the hard functions H^. The corresponding one-loop expressions have been given in (25). At 
that scale, they are free of large logarithms and so can be reliably computed using perturbation 
theory. We then evolve the hard functions down to an intermediate hard-collinear scale /ij ~ 
iy/m b AQ CD by solving the renormalization-group equation 

H i:j (n -p,//) = 27j(n -p,//) H^n -p,//) , (35) 



dln/i 



where 7j is the anomalous dimension of the semileptonic current in SCET. 

ii) Next, we start from a model for the shape function S(u,fio) at some low scale /io = 
few x Aq CD large enough to trust perturbation theory. Such a model could be provided by 
a QCD-inspired approach such as QCD sum rules or lattice QCD, or it could be tuned to 
experimental data. We then solve the integro-differential evolution equation 



d 



dln/i 



S(cu,n) = —J duo'^si 1 ^-, fi) S(u', n) (36) 



to obtain the shape function at the intermediate scale /ij. 

iii) Finally, at the scale /ij we combine the results for the hard functions and for the shape 
function with the jet function J in (26), which at that scale is free of large logarithms and 
so has a reliable perturbative expansion. The dependence on the matching scales Hh and /z; 
cancels in the final result (to the order at which we are working). 

We now discuss these three steps in detail. 
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4.1 Evolution of the hard functions 



At one-loop order, the anomalous dimension 7j for the SCET current is twice the coefficient 
of the 1/e pole in the renormalization factor Zj in (24). More generally [19, 29], 



7j(n • p, u) = -r cusp (a s ) In + 7' (<x, N 



n ■ p 



IT 



-In 



n ■ p 



+ ... , 



(37) 



where T 



cusp 



Cf^s/^ + • • • is the universal cusp anomalous dimension governing the ultra- 
violet singularities of Wilson lines with light-like segments [30]. The exact solution to the 
evolution equation (35) can be written as 



Hij{n ■ p, in) = Hijin ■ p, fj, h ) exp U H (n ■ p, fi h , , 



where 



U H (n ■ p, n h , Hi) = 2 



da 
P{a) 



Tcusp (c) 



In 



n ■ p 
Hh 



da' 



+ i(a) 



(38) 



(39) 



and f3(a s ) = da s /d\rtn is the QCD (5- function. Defining as usual 



Tcusp (^s 



emit 1 

n=0 



f3(a s ) = -2a s J2Pn(^) 



(40) 



n=0 



and similarly for all other anomalous dimensions, setting r\ = a s (/ii) / 'a s (fJ>h) > 1> an d ex- 
panding the evolution function to 0(a s ), we obtain 



e u H (n-p,n h ,M) _ e v H (fi h ,m) ( n ' P \ 130 



Hh 



A In ri 



a s (fih) T /r\ PA fi-p 

4tt Aj Vr Po) 1 1 ; Hh . 



, (41) 



where 



V H (Hh,Hi) 



— — 1 lnri + — In n - — - — ) (ri - 1 - lnrj 

V r i / 2 Po \ r o Po/ 



^ In n + 

Po 



(ri - 1) a s {n h ) 



(42) 



Let us briefly explain the structure of resummed perturbation theory in applications with 
Sudakov double logarithms. In renormalization-group improved perturbation theory the pa- 
rameter r\ is treated as a quantity of 0(1). The terms proportional to l/a s {iXh) in Vh resum 
the leading, double logarithmic terms to all orders in perturbation theory. The remaining 
0(1) terms in Vh contribute at leading, single-logarithmic order. Note that these effects are 
not suppressed by any small parameter. We therefore refer to the combination of these two 
terms as the "leading order". At next-to-leading order, the corrections proportional to the 
coupling a s (fih) are included. In our case, the only piece missing for a complete resumma- 
tion at next-to-leading order is the 0(a s ) contribution to Vh, which is independent of the 
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kinematic variable n ■ p and vanishes for ^ — > fi h . To compute these terms would require to 
calculate the cusp anomalous dimension to three loops and the anomalous dimension 7' to two 
loops. The fact that the corresponding expansion coefficients are unknown implies a universal, 
process-independent small uncertainty in the normalization of inclusive 5-decay spectra in the 
shape-function region. We stress, however, that this uncertainty cancels in all ratios of decay 
distributions, even between B — > X u l~u and B — > X s j spectra. 

In phenomenological applications of our results the intermediate scale ^ will be of order 
m c , because a restriction to hadronic invariant masses below the charm threshold is used to 
separate B — > X u l~v from B — > X c l~v decays. It is then appropriate to perform the running 
between \i h and /ij in a theory with rif = 4 light quark flavors. The relevant expansion 
coefficients are T = f,T 1 = ^- fir 2 , % = -f , and (3 = f , ft = ±fk 



4.2 Evolution of the shape function 



At one-loop order, the anomalous dimension for the shape function is twice the coefficient of 
the 1/e pole in the renormalization factor Z s . From (31), we obtain 



C F a s 



7T 



2 In 



A- uj 



1 ) 8(u - uj') - 2 



8(u' 



oj 



UJ' 



UJ 



(43) 



Let us briefly compare this result with previous calculations of this anomalous dimension pub- 
lished in the literature. In [16], Balzereit et al. obtained a result that almost agrees with ours, 
except for the factors of 2 in front of the logarithmic term and of the plus distribution. (Their 
expression is, however, given in a rather different form.) Aglietti and Ricciardi tried to extract 
the renormalized shape function by matching QCD directly onto HQET, without including a 
jet function [31]. They obtained results for the renormalization factor Zs (and hence for the 
anomalous dimension) in two different regularization schemes. Their expression obtained in 
dimensional regularization disagrees with our findings. Bauer et al. computed the anomalous 
dimension of the shape function using SCET [10]. They interpreted the renormalization factor 
in (31) as a distribution in uj rather than u/ and studied the renormalization-group equation 
for the short- distance coefficients Hij ■ J instead of that for the shape function. If this is done, 
the logarithm in the anomalous dimension in (43) must be replaced by ln[/x/ (Ahad + ^-OL where 
Ahad is the lower cutoff on the integral over ui (see the discussion in Section 3.3). The value of 
this cutoff is process dependent and set by kinematics or by virtue of some experimental cut. 
In the shape-function region, A ha d = 0(Aq CD ). Instead, in [10] the combination (A ha d +oj') is 
identified with the 6-quark mass, which is incorrect. 

The evolution equation (36) can be solved analytically using a general method developed 
in [15]. It is convenient to change variables from uj to uj = A — uj G [0,oo[ and denote 
S(uj) = S(A — uj). The renormalization-group equation then reads 



d 



duj'^ s (uj,uj',/j,) S(u',fi) , 



din// j 
where the anomalous dimension can be written in the general form 



r cU s P (« s ) ln^ + l(a s ) 

UJ 



S(uj — uj') + 2Q(uj, uj', a s 



(44) 



(45) 
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The logarithmic term containing the cusp anomalous dimension has a geometric origin. Since 
the heavy-quark field h(x) in HQET can be represented as a Wilson line along the v direc- 
tion, the field "K(x) entering the SCET formalism contains the product of a light-like Wilson 
line (along n) and a time-like Wilson line (along v), which form a cusp at point x. The 
shape function contains two such cusps. According to the renormalization theory of Wilson 
lines with light-like segments, each cusp produces a contribution to the anomalous dimension 
proportional to r cusp ln/z [30]. The one-loop coefficients of the remaining terms in (45) are 



7o = -2C, 



B{u - Co 1 ) 



Co — Co' 

The general solution of (44) can be obtained using the fact that on dimensional grounds 



(46) 



roc 

I duo' g(Cj,Qj\a s )(Cj') a = Cu a ^(a,a s ) , 
Jo 



(47) 



where the function 5F only depends on the exponent a and the coupling constant. We set 
jF(0, a s ) = by definition, thereby determining the split between the terms with 7 and Q in 
(45). The integral on the left-hand side is convergent as long as Re a > — 1. At one-loop order 
we find from (46) 

ex r 

?(a,a s ) = T ^ [i>(l + a)+ lE \ +••• , (48) 

where ip(z) is the logarithmic derivative of the Euler T function. Relation (47) implies that 
the ansatz [15] 

/ £ \ r+2g(ji,no) 

exp£7 5 (T,/z,/! ) (49) 



f(u,H,Ho,T) = — 



with 



(50) 



Us(t,/i,iiq) = -2 



da 



g(n, Ha) + l(a) + 3{t + 2g(p a , /i ), a) 



provides a solution to the evolution equation (44) with initial condition /(a), // , Ho, T ) — 
[uj/ HoY at some scale /i . Here \i a is defined such that a s (/i a ) = a, and r can be an arbitrary 
complex parameter. Note that g(/i, /x ) > if /i > Ho- We now assume that the shape function 
S(u,fio) is given at the low scale /i and define its Fourier transform with respect to ln(d>/// ) 
through 

1 f°° ( uj\ U 

S(cu,Ho) = 7T / dtSo(t)i — ) . 

The exact result for the shape function at a different scale jj, is then given by 



(51) 



^,/i) = - 



dtS (t) f(u;,n,n ,it) 



(52) 
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With the help of this formula, it is straightforward to derive explicit expressions for the 
evolution of the shape function from the hadronic scale liq up to the intermediate scale /ij at any 
order in renormalization-group improved perturbation theory. Setting r 2 = a s (fio)/a s (fii) > 1, 
we obtain for the evolution function at leading order 



f(u,Lli,Ll ,it) 



e Vs(w>«>) 

V/4> 



r(i +it) 



r(l + ft + glnr 2 )' 



(53) 



where 



Vs(/Xi,/x ) = 



47T 



a s (/xo) 



r 2 



In r 2 

r 2 



lnr 2 - ^ lnr 2 + (r 2 - 1) a s (/i) 
Po Po 



(54) 



This result is valid as long as (r /po) l nr 2 < 1) which is the case for all reasonable parameter 
values. Missing for a resummation at next-to- leading order are the 0(a s ) contributions to Vg, 
which vanish for ii — > fj, . Since these corrections have an unknown dependence on t via the 
two-loop contribution to the function 9 r (o; s ,a), they will affect the Co dependence of the final 
result. There are also some known 0(a s ) corrections to (53) proportional to m(a>//x ), which 
we have omitted for consistency. For all practical purposes, given the intrinsic uncertainties in 
our knowledge of the shape function, it will be sufficient to use the equations given above. As 
mentioned earlier, we typically have \ii ~ m c , and so the running between Lii and li should 
be performed in a theory with rif — 3 light quark flavors. The relevant expansion coefficients 
are then r = f,F 1 = ^f- f vr 2 , 7o = -|, and (3 = 9, (3 X = 64. 

The leading-order result presented above can be simplified further. When (53) is inserted 
into (52), the integration over t can be performed analytically. Setting 77 = (T /f3 ) lnr 2 > 0, 
the relevant integral is 



1 r°° 

1 = yJ 



T(l + it) 

r(i + it + r]) ' 



(55) 



So(t) 



f 

Jo 



(56) 



where 

' duj' *. , 1 uj 

is the Fourier transform of the shape function as defined in (51). The integrand of the t- 
integral has poles on the positive imaginary axis located at t — in with n > 1 an integer. For 
Co < Co 1 the integration contour can be closed in the lower half-plane avoiding all poles, hence 
yielding zero. For Cj > Co' we use the theorem of residues to obtain 



1= dCo' R(Co,Co') S(Co' , liq) , 
Jo 



where 



1 J uoj^K Co) T{] + 1 



(57) 



(58) 
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Note that R{Cj,Cj') — > 5{Cj — Cj') in the limit rj — > 0, corresponding to /x^ — * /i , as it should 
be. Our final result for the shape function at the intermediate hard-collinear scale, valid at 
leading order in renormalization-group improved perturbation theory, can now be written in 
the simple form (valid for /Xj > /x , so that r\ > 0) 

S(Cj, = e^o) 1 W' 1 ?) , (59) 

rfa) Jo n - w') 1 ^ 

with as given in (54). A similar analytic result for the renormalized shape function was 
obtained in [16] using a different strategy to solve the evolution equation for the shape function. 
However, these authors miss a factor 2 in the expression for 77, and we disagree with their 
expression for the function Vs. 

From the above equation one can derive scaling relations for the asymptotic behavior of the 
shape function for Cj — > and Cj — > 00 (corresponding to uj — > A and — > —00). If the function 
S(Cj, fxo) at the low scale /xo vanishes proportional to u)^ near the endpoint, the shape function 
at a higher scale /Xj > /x vanishes faster, proportional to Cj^ +ri . Similarly, if S(u, /xo) falls off like 
Cj~^ for Cj — > 00, the shape function renormalized at a higher scale vanishes like ct)~ mm ( 1 '^ +?? . 
Irrespective of the initial behavior of the shape function, evolution effects generate a radiative 
tail that falls off slower than 1/oj. This fact implies that the normalization integral of S(u),fx) 
as well as all positive moments are ultra-violet divergent. The field-theoretic reason is that the 
bilocal shape-function operator in (9) contains ultra-violet singularities as z_ — > 0, which are 
not subtracted in the renormalization of the shape function. The situation is analogous to the 
case of the B- meson light-cone distribution amplitude discussed in [15, 32]. These divergences 
are never an obstacle in practice. Convolution integrals with the shape function are always 
cut off at some finite value of Cj by virtue of phase-space or some experimental cut. 



5 Properties of the shape function 

In this section we discuss how moments of the shape function are related with HQET param- 
eters. This will lead us to propose a new, physical scheme for defining a running heavy-quark 
mass, which is most appropriate for the study of inclusive spectra in the shape-function region. 
We will also present a model-independent result for the asymptotic behavior of the renormal- 
ized shape function (defined in the MS scheme), finding that it is not positive definite. 

Most of our discussion in this section is phrased in terms of the original (unhatted) shape 
function S(u,/j,). At the end, we formulate the resulting constraints on the function S(u,fj). 



5.1 Shape-function moments in the pole scheme 

Naively, ignoring renormalization effects, the moments Mn = du oj n S(oj) are given by 
hadronic parameters defined in terms of 5-meson matrix elements of local HQET operators 
[1]. In particular, M = 1 fixes the normalization of the shape function, Mi = vanishes 
by the HQET equation of motion, and M 2 = — Ai/3 is determined by the matrix element of 
the kinetic-energy operator. The vanishing of the first moment is connected with the implicit 
definition of the heavy-quark pole mass built into the HQET Lagrangian via the equation 
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of motion iv ■ D h = 0. These moment constraints have been implemented in various model 
parameterizations for the shape function suggested in the literature [3, 33, 34]. Typically, one 
makes an ansatz for the shape function depending on a few HQET parameters such as A and 
Ai, and determines the values of these parameters from a fit to experimental data. 

Beyond tree level, all moments M N with N > receive ultra-violet divergences from the 
region uj — > — oo (or Cj — > oo). However, as we have mentioned above, the values of ui needed 
for the description of physical decay rates are always restricted to a finite interval. It is thus 
sufficient for all purposes to define the moments of the renormalized shape function as 

,A 

Mjv(A uv ,/i)= / dujuj N S{u,ii). (60) 

</-A uv 

The dependence of these moments on the renormalization scale fi is controlled by the evolution 
equation (36). In addition, the moments depend on the lower cutoff on the u integral. The 
choice of Auv is a matter of convenience, and so we are free to pick a value that is numerically 
(if not parametrically) large compared with Aq CD . In this case, as we will now show, the 
dependence on Auv can also be controlled using short-distance methods. 

For sufficiently large values of Auv if is possible to expand the moments M^(A\jy, n) in 
a series of 5-meson matrix elements of local HQET operators. If for simplicity we set T = 1 
in the shape-function operator (which is legitimate, since the Dirac structure is unaltered in 
HQET), the operators in question are Lorentz-scalar, "leading-twist" operators containing 
h...h [1, 2]. These are the operators that mix with h(in ■ D) N h under renormalization. 
It is straightforward to find the corresponding operators of a given dimension. The unique 
dimension-3 operator is hh. The two operators of dimension 4 are hin ■ Dh (class-1) and 
hiv ■ Dh (class-2). The class-2 operator vanishes by the HQET equation of motion. The 
possible dimension-5 operators are 

class-1: h(in-D) 2 h, h(iD ± ) 2 h, 

(61) 

class-2: h(iv-D) 2 h, hin-Div-Dh, hiv- Din- Dh, 

where again the class-2 operators vanish by the equation of motion. Moreover, it follows from 
the Feynman rules of HQET that the two class-1 operators do not mix under renormalization, 
so the operator h(iD±) 2 h can be ignored. From dimension 6 on the situation is more com- 
plicated, because several class-1 operators exist that can mix with h (in ■ D) N h. For N = 3 
these are of the form hiD Gh or J2 q h . . . qq . . . h, where we omit Lorentz and color indices. 
We will restrict our discussion to operators of dimension less than 6. 

For the operator product expansion of the moments in (60) we need the forward matrix 
elements 

of the leading-twist operators between 5-meson states in HQET. Using the equation of motion, 
it can be shown that (hh) = 1, (hin ■ Dh) =0, and (h(in ■ D) 2 h) = — Ai/3 [1]. We can thus 
write an expansion of the form 

3 ■ 



M N (A vy ,fi) = A^ v { K { N) (A vy ,f,) + K ( 2 N } (Auv,/i) • + O 



Aqcd 



A 



uv 



(63) 
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This expansion is useful as long as the cutoff Auv is chosen much larger than the typical 
hadronic scale characterizing the matrix elements of the local operators. The matching co- 
efficients in this relation can be calculated using on-shell external 6-quark states with 
residual momentum k. For operators of dimension up to 5 it suffices to calculate two-point 
functions. (Three and four-point functions would have to be considered at dimension 6.) We 
first evaluate the moments of the renormalized shape function in (33), finding at one-loop 
order 



Mr ton (Auv,/i) = (n-A;) 



N 



C F a s 



7T 



In 



2 Auv + n ■ k i , Auv + n ■ k i n 2 



+ ln 



// 



+ 24 



N 



N 



C'i <Y. S 1 j 1 + 2 ln A uv + n ■ k ^ 2 



7T 



J 



A 



uv 



n ■ k 



1 



(64) 

i //•/,' / i 

j=l J \ ^ 1=3 

We then expand this result in powers of n ■ fc/A uv . Keeping the first three terms in the 
expansion, we obtain 



M parton (A uv ,/i) = l 



C F a s 



7Y 



In' 



A 



uv 



C F a s 



7Y 



Mr ton (Avv^) = n-k 



1 - 



n ■ k 
A^v" 

C F a, 



21n 
In 



+ ln 


Auv 


7T 2 \ 


n 




Auv 


*■) 


(n ■ kf 


n 


+ A 2 



7T 



2 ^-UV 



fi 2 



-In 



A 



uv 



7T 

24" 1 



M 2 parton (A uv ,/i) 



_ C F a s 

7T 

(n ■ kf 

7T 



A 



uv 



_ 21n ^ + 1 U^21n^ + 



7T 



In 



2 -^-UV 



A 



uv 



// 



(65) 



// 24 2 



A 2 
A uv 



ln^ • /,•/,• A, v (-■>\n ±± • 3H . 



In the next step, we calculate the one-loop matrix elements of the local operators h (in- D) N h 
between heavy-quark states with residual momentum k. The relevant diagrams are the same 
as in Figure 2, where now the black square represents the local operators. Keeping v ■ k 
non-zero to regularize infra-red singularities, we obtain for the bare matrix elements 



(h (in ■ D fh) = (n-k) 



v 1 I - (-2v ■ k) 2 < 



N 

E 



N 
3 



'2v ■ k\ j 



(An) 2 -* 

x(j-l-e)T(j-e)T(2e-j)}. 



n ■ k 



(66) 



While the individual diagrams are infra-red divergent, taking the limit v ■ k — > in the 
sum of all contributions is possible without encountering singularities. Then the one-loop 



19 



contributions vanish, and the matrix elements simply reduce to their tree-level values. In 
other words, the one-loop contributions correspond to a mixing with class-2 operators, whose 
hadronic matrix elements vanish by the equations of motions. It follows that in (65) we must 
identify (n ■ k) n — > (h(in ■ D) n h). Substituting the results for the HQET matrix elements 
given earlier, we obtain for the Wilson coefficients of the first three moments 

K io) = 1 _C £ a 1 ( 2 Auv Aw f\ ( 0) = / Anv _ 1 

/fW = ^f21n^-iy K? = -2^\n^ : (67) 



7t \ H / 7T [I 

= - l n , K f = i _ 9^1 ( ln 2 _ 2 In + ± _ I 

° 7T // ' 2 7T V ^ /i 24 2 

At tree level, this reproduces the naive moment relations mentioned at the beginning of this 
section. Beyond tree level, the moments get corrected by calculable short- distance effects, 
which can be controlled using fixed-order perturbation theory as long as the ratio Auv/ A* is °f 
0(1). In particular, the renormalized first moment no longer vanishes, but is proportional to 
the cutoff Auv up to small power corrections. 

As mentioned earlier, the value of the first moment is connected with the definition of the 
heavy-quark mass (see also [26, 28]). The first moment of the renormalized shape function can 
be made to vanish to all orders in perturbation theory by choosing an appropriate scheme for 
the definition of m b . So far our calculations have assumed the definition of the heavy-quark 
mass as a pole mass, ml° le , which is implied by the HQET equation of motion iv ■ D h = 0. 
Results such as (67) are valid in this particular scheme. A more general choice is to allow 
for a residual mass term 5m in HQET, such that iv ■ D h = 5m h with 5m = O(Aqcd) [35]. 
It is well known that the pole mass is an ill-defined concept, which suffers from infra-red 
renormalon ambiguities [36, 37]. The parameter A po i e = Mb — ml° le , which determines the 
support of the shape function in the pole-mass scheme, inherits the same ambiguities. It 
is therefore advantageous to eliminate the pole mass in favor of some short-distance mass. 
For the analysis of inclusive _B-meson decays, a proper choice is to use a so-called low-scale 
subtracted heavy-quark mass m^iif) [38], which is obtained from the pole mass by removing 
a long-distance contribution proportional to a subtraction scale [if = few x Aq CD , 



m b 



= m b (fif) + fif g(a s (fj,),^j = m b (/i f ) + 5m. (68) 
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As long as m b (nf) is defined in a physical way, the resulting perturbative expressions after 
elimination of the pole mass are well-behaved and not plagued by renormalon ambiguities. 
Replacing the pole mass by the physical mass shifts the values of n ■ k and uo by an amount 
5m, since n ■ (ml° le v + k) — m b (fif) + (n ■ k + 5m), and because the covariant derivative in 
the definition of the shape function in (13) must be replaced by in ■ D — 5m [35]. At the 
same time, A po i c = A(/i/) — 5m, where A-(uf) = Mb — m b (uf) is a physical parameter. Note 
that this leaves the parameter Cj = A — uj and hence the shape function S{Cj,u) invariant. 
This follows since A pole — oj vo \ c = A(fij) — (u po \ c + 5m), where u; po i e denotes the value in the 
pole-mass scheme used so far. 
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5.2 Shape-function moments in a physical scheme 

From now on we will adopt a mass scheme defined by some specific choice of 5m. Let us denote 
by uj = cjp i c + 5m the value of the light-cone momentum variable in that scheme and define 
"physical" moments M^ hys as in (60), but with all parameters replaced by their values in the 
new scheme, in particular A = A(fif). Then the expressions for the moments in (64) and (65) 
change according to the replacements n-k — > n-k + 5m and Auv — > Auv — 5m everywhere. We 
now choose 5m such that the first moment vanishes, thereby defining a low-scale subtracted 
heavy-quark mass (with /if — Auv) to all orders in perturbation theory We will refer to this 
mass as the "shape-function mass" mf F . This is a "physical", short-distance mass in the sense 
that it is free of renormalon ambiguities. (However, the definition of the shape-function mass 
depends on the renormalization scheme used to define the shape function.) From (65) and 
(68), it follows that at one-loop order 
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Note that after introduction of the shape-function mass the coefficients K n in the operator- 
product expansion for the moments in (63) vanish by definition. However, to first order in a s 
the values for the coefficients K$ and K^f 1 of the zeroth and second moments given in (67) 
remain unchanged, since 5m = 0(a s ). This would no longer be true for the coefficients of 
higher moments. 

The shape-function mass can be related to any other short-distance mass using pertur- 
bation theory. For instance, at one-loop order its relations to the potential-subtracted mass 
introduced in [39] and to the kinetic mass defined in [40, 41] read 
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(70) 



Note that, in addition to the dependence on the subtraction scale ///, the shape-function mass 
depends on the scale /i at which the shape function is renormalized. While it is natural to set 
fi = /if, as we did here, this is not necessary. Given a value for the shape-function mass for 
some choice of scales, we can solve (69) to obtain its value for any other choice, using the fact 
that the pole mass is scale independent. 

Proceeding in an analogous way, we can use the second moment to define a physical 
kinetic-energy parameter, commonly called This quantity can be used to replace the 
HQET parameter Ai, which like the pole mass suffers from infra-red renormalon ambiguities 
[42]. At one- loop order, we obtain 
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Taking the ratio of M2 hys and MQ hys has the advantage of eliminating the double logarithmic 
radiative corrections from this expression. Our definition is similar to the running parameter 
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fi\ defined in the kinetic scheme [40, 41]. At one-loop order, the two parameters are related 
by 

C F a s (fi f ) 
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(72) 



Given a value for the kinetic energy in the shape-function scheme for some choice of scales, 
we can solve (71) to obtain its value for any other choice, using that Ai is scale independent. 

Similarly, each new moment of the renormalized shape function can be used to define a 
new physical, scale-dependent parameter A^^Auy, jj) = M^ ys (A\jy , fi) / MQ hys (Auy , fi) , which 
coincides with the corresponding HQET parameter = (h(in ■ D) N h) at tree level, and 
which beyond tree level is related to HQET parameters through well-controlled perturbative 
expressions. Obviously, the presence of power divergences implies that higher moments are 
progressively less sensitive to HQET parameters, since they are dominated by the perturbative 
terms of order a s A^ Y . 



5.3 Moments of the scheme-independent function S(u>,fi) 

It will be useful to rewrite the moment relations derived above in terms of the variable Co = 
A — uj, which is invariant under redefinitions of the heavy-quark mass. Defining a new set of 
scheme-independent moments 
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we obtain 
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where the parameters A(/i/,/i) = Mb — mf F ([if,fi) and fi^(fif,fi) should be considered as 
known physical quantities. Using the relations in the previous subsection, we have 
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where /x* denotes the scale at which initial values for the two parameters are obtained, for 
instance using relations such as (70) and (72). These relations are particularly simple if one 
chooses fif = fi. 
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In the relations above we have eliminated the unphysical HQET parameter Ai in favor of 
the physical parameter jj^ defined in the shape-function scheme. At first sight, this seems 
to threaten the convergence of the operator product expansion. For instance, the term pro- 
portional to fx\ in the expression for the zeroth moment Mo in (74) contains a leading-power 
perturbative contribution of order a^(/i), and similar contributions would arise from all other 
terms in the expansion. These contributions would have to be subtracted from the Wilson 
coefficient of the first term, if this coefficient were computed to two-loop order. The overall 
convergence of the operator product expansion is unaffected by this reorganization of pertur- 
bative corrections. 



5.4 Asymptotic behavior of the shape function 



The fact that for sufficiently large values of the cutoff the moments of the shape function 
can be calculated using an operator-product expansion implies that a similar expansion can 
be used to obtain a model-independent description of the asymptotic behavior of the shape 
function. Taking the derivative of the zeroth moment M in (73) with respect to /if, one 
obtains 
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This relation can be trusted as long as /if ^> Aq CD . It allows us to determine the behavior of 
the shape function for large values of a). From (74) we find at one-loop order 
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(77) 



The precise definitions of A and [i\ are not specified at this order. (Note that the shape 
function cannot depend on the value of the cutoff ///.) We have checked that this asymptotic 
behavior of the shape function is consistent with the evolution equation (59) when expanded 
to first order in a s . 

Relation (77) is a model-independent result as long as Co 3> A QCD . We stress the remark- 
able fact that this radiative tail of the shape function is negative, in contrast with the naive 
expectation based on a probabilistic interpretation of the shape function as a momentum dis- 
tribution function. The point is that the definition of the renormalized shape function requires 
scheme-dependent ultra-violet subtractions. From (77) it follows that the shape function must 
have a zero, which for sufficiently large \i is located at a value ujq ~ A + ji/ y/e. 



6 Recapitulation 

We have now completed the conceptual part of this paper. Before turning to phenomenological 
applications, let us briefly summarize the discussion so far. The hadronic physics governing the 
inclusive semileptonic decay B — > X u l~v is encoded in the structure functions Wi appearing in 
the Lorentz decomposition of the hadronic tensor W^ u in (17). In the shape-function region, 
only two combinations of these functions are required at leading order in Aqcd/^6- They 
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follow from the factorization formula (14) using the explicit results for the hard functions 
and the jet function J derived in this paper. Explicitly, we obtain at next-to-leading order in 
renormalization-group improved perturbation theory 
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where the dots represent power corrections in Aqcd/^6- In these expressions y — n • p/rrib 
is a partonic scaling variable, while = m^n • Pr — CS) depends on the hadronic variable 
n • Pr = Eh — \Ph\- In ah our results, nib denotes the heavy-quark mass defined in the 
"shape-function scheme" introduced in Section 5 (see also Section 9 below). The jet function 
J(ftu,y, Hi) at an intermediate hard-collinear scale \ii ~ ^JTf^K^p^ can be calculated in fixed- 
order perturbation theory. The relevant expression valid at one- loop order is given in (27). The 
above results contain a variety of renormalization-group functions, which arise in the solution 
of evolution equations discussed in Section 4. The explicit form of the Sudakov exponent V H 
can be found in (42). This function is independent of the kinematic variables y and jjfc. In 
addition, we need 
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For simplicity, we have identified the high-energy matching scale /ih introduced in Section 4 
with the heavy-quark mass m b . Our results are formally independent of the precise choice of 
Hh ~ m b . The numerical effect of the residual \ih dependence remaining after truncation of 
the perturbative expansion has been studied in [29] and was found to be small. 

The function S(u, fii) in (78) is the shape function after the transformation of variables 
from 00 to Co = A — uo. The limits of integration for the variable Co (i.e., < Co < n ■ Ph) are 
set by hadronic kinematics and are independent of the definition of the heavy-quark mass. 
The shape function is a non-perturbative object, which at present cannot be predicted from 
first principles. It enters our results (78) renormalized at the intermediate hard-collinear scale 
Hi. In (59), we have presented an analytic formula (valid at leading order in renormalization- 
group improved perturbation theory) that relates the shape function at a high scale to the 
shape function renormalized at a low hadronic scale. Many properties of the shape function 
that were so far unknown have been derived in Section 5. In particular, we have given explicit 
formulae relating the moments of the shape function to HQET parameters, and we have proved 
that the shape function has a negative tail for large values of Co, whose explicit form can be 
calculated using an operator product expansion. These new insights about the shape function 
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will be very helpful in constructing a realistic model for the function S(u,fj,i), which can then 
be refined by tuning it to experimental data such as the photon energy spectrum in inclusive 
B — > X s 7 decays. 



7 Differential decay rates and spectra 



In the shape-function region considered in this work, the hadronic tensor is most naturally 
expressed in terms of the variables n • Pr and n ■ p, where p = m b v — q = Pu — At> is the 
momentum of the final-state hadronic jet in the parton picture. This would remain true if 
we worked to higher order in the collinear expansion. It is thus useful to derive expressions 
for the decay rates in terms of these variables. Our theoretical results are valid as long as 
n ■ P H can be considered as being of order a hadronic scale (say, a few x Aq CD ), whereas 
n ■ p is integrated over a domain of order m b ^> Aq CD . It is this integration which provides a 
sampling over sufficiently many hadronic final states needed to ensure quark-hadron duality. 
The duality hypothesis underlies any description of inclusive decay rates using short- distance 
methods [43]. 

Under these conditions, it is appropriate to describe the distribution in n ■ p in terms of 
a partonic scaling variable y = n ■ p/m b , while the distribution in the orthogonal light-cone 

— * 

component is described in terms of the dimensionful hadronic variable P + = n-Pn = Eh — \Ph\- 
At leading order in Aqcd/^6, we obtain from [18] the triple differential decay rate 
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where x — 1 — x, and x = 2Ei/m b is a scaling variable proportional to the energy of the 
charged lepton measured in the S-meson rest frame. r tree = G 2 F \V ub \ 2 (ml° lc ) 5 / (192n 3 ) denotes 
the tree-level expression for the total B — > X u l~v decay rate obtained at leading order in the 
heavy-quark expansion. Phase space is such that 



< P+ < M B - 2E t = m b x + A . 



P + -A 
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(81) 



The structure functions Wi have support only for y > 0, see (27). If the lepton scaling variable 
x is integrated over a domain of order unity (meaning that Ei is integrated over a domain of 
order m b ^> Aqcd), one can replace the second condition by0<x<y<lat leading power 
in Aqcd/"^?)- If, on the other hand, the lepton energy is restricted to be close to its kinematic 
limit, Ei PS Mb/2, then x = 0(Aqcd/"^6), and at leading order the rate (80) can be simplified 
to 
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with < y < 1 and < P + < M B - 2E X . 

It will be useful to develop some intuition for the light-cone momentum variables. In 
general, the hadronic tensor can be described in terms of the quantities 
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Figure 3: Hadronic phase space for the light-cone variables P_ and P + (left), 
and theory phase space for m& = 4.8 GeV (right). The scatter points indicate 
the distribution of events as predicted by the model of [18]. In each plot the 
solid line separates the regions where sh < Mp. (dark gray) and sh > Mp. 
(light gray), whereas the dashed line corresponds to P + = M^/Mb- The 
dotted line in the first plot shows the contour where q 2 = (Mb — Mb) 2 - 

whose true phase-space is M 2 /P_ < P+ < P- < Mb, corresponding to a triangular region in 
the (P-, P + ) plane with a tiny portion near the P_ axis left unpopulated. The variable P_ is 
related to our parton variables by P_ = n ■ p + A = m^y + A. In our theoretical description 
based on quark-hadron duality P + starts from 0, while the small region with P_ < A is 
left unpopulated. This is illustrated in Figure 3. Contours of constant hadronic or leptonic 
invariant mass in the (P_, P + ) plane are easy to visualize, since 

s H = P 2 =P + P_, q 2 = (M B - P + )(M B - P-) (84) 

are given by simple expressions. The solid and dotted lines in the left-hand plot in Figure 3 
show the contours where sh = Mp. and q 2 = (Mb — Mb) 2 , respectively, which can be used 
to separate B — > X u l~v events from semileptonic decays with charm hadrons in the final 
state. The dashed horizontal line shows the maximum allowed value of P + when a cut Ei > 
(M 2 B — M d )/(2M b ) is applied to the charged-lepton energy, which implies P + < M^/M B . 
This cut is another way of eliminating the charm background. In the right-hand plot, we 
indicate the density of events in theory phase space obtained using the model of [18]. 1 The 

1 While not rigorously implementing shape-function effects beyond tree level, the model of [18] has the 
advantage that it interpolates between the shape-function region and the remainder of phase space, where 
a local operator product expansion can be employed. On the contrary, our more rigorous discussion here is 
limited to the region of hard-collinear jet momenta. We believe that the scatter plot shown in the figure 
provides a reasonably realistic impression of the population in phase space. 
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vast majority of events is located in the shape- function region of small P + and large P_. 

In the remainder of this section, we present exact analytic results, valid at leading order in 
Aqcd/^?) and at next-to- leading order in renormalization-group improved perturbation theory, 
for a variety of spectra in B — > X u l~u decays. They are obtained by using our results (78) 
in conjunction with the expressions for the differential rates in (80) or (82), as appropriate. 
Our strategy will always be to integrate over the scaling variable y before integrating over the 
hadronic variable P + , changing variables from P + to p 2 = mf ) (P + — Co). In this step, one must 
carefully evaluate the effect of the star distributions contained in the jet function J(p 2 , y, Hi), 
using the definitions in (19). The integral over the shape-function variable Co is left until the 
end, so that our formulae can be evaluated once an explicit form for the shape function is 
assumed. We will always present fractional decay rates normalized to the total inclusive rate 
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where the dots represent higher-order perturbative corrections as well as power corrections 
of order (AqcdAt"^) 2 and higher. This procedure offers the advantage of eliminating the 
strong sensitivity to the heavy-quark (pole) mass. Replacing r trec by r tot adds a contribution 
(27T 2 — y) to the coefficient of the hard correction to the function W\ in (78). Our predictions 
for normalized rate fractions can be turned into predictions for absolute rates with the help 
of an independent theoretical prediction for the total B — > X u l~u rate obtained using a local 
operator product expansion. 

The integrals over the parton variable y encountered in our analysis can be reduced to a 
set of master integrals defined as 
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where 6 > — 1 and z < 1 are a arbitrary real numbers. The integral J 4 can be expressed in 
terms of the incomplete (3 function B(z, a, 6) and the Lerch transcendent $(2, a, 6) as 



h(b, z)=lnz B(z, 1 + 6, 0) - z 1+h 2,1 + 6). 
We note the useful relations 

h{\ + b,z) = I 4 (b, z) - J 2 (6, z) , J 4 (6, 1) = -^'(1 + 6) , 
where ip'(z) is the derivative of the Euler ip function. 
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7.1 Charged-lepton energy spectrum 



As a first application, we study the distribution of the charged-lepton energy near the kine- 
matic endpoint. Specifically, we assume that M B — 2Ei is of order a hadronic scale. Starting 
from the triple differential rate in (82), we obtain for the normalized energy spectrum 
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Here 
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Using the results for the master integrals in (86), one readily derives the explicit formulae 
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At leading power in Aqcd/^6 the heavy-quark mass in the denominator of the prefactor on 
the right-hand side of (89) can be replaced by m b + uj = M b — Cj } which removes any sensitivity 
to the definition of m b . This replacement can indeed be justified by studying power corrections 
to the shape function [44, 45]. 

All our results for decay rates will have a similar structure, but the definitions of the 
functions T, /„, and 77 will be different in each case. The product T e VtI resums the leading 
double and single-logarithmic corrections to all orders in perturbation theory. The tree-level 
result can be recovered by setting Vh — and a = 0, in which case T(0) = 1, and the spectrum 
is simply given in terms of an integral over the shape function [1]. The next-to- leading order 
terms can be divided into an ^-independent hard function 77(a), whose structure follows from 
the form of the hard corrections in (78), and a sum of hard-collinear radiative corrections, 
which follow from the integration over the jet function. 



28 



Using the above result, it is straightforward to calculate the fraction F E = T(Ei > E )/T tot 
of all B — > X u l~u events with charged-lepton energy above a threshold E . Defining A E = 
M B - 2E , we find 
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The fraction F^A^) is given in terms of a weighted integral over the shape function, with a 
weight factor of order Aqcd/^6 that vanishes at the upper end of integration. As a result, 
only a small fraction of events is contained in the lepton endpoint region. 



7.2 Hadronic P + spectrum 



A cut on the charged-lepton energy restricts the variable P + to be less than A E . On the 
contrary, however, a cut on P + does not restrict the lepton energy to be in the endpoint 
region. In fact, according to (81) the integration over x can be taken to run from to 1 at 
leading power if P + is small. It follows that the fraction of events with P + < A E samples 
the same hadronic phase space as the lepton-endpoint cut, but it contains significantly more 
events. Such a cut therefore offers an excellent opportunity to determine the CKM matrix 
element \V ub \. 

The calculation of the fraction Fp(Ap) = T(P + < Ap)/T tot of all B — > X u l~v events with 
hadronic light-cone momentum P + below a threshold Ap starts from the expression for the 
triple differential rate in (80). We integrate over x and y in the range < x < y < 1 before 
integrating over P + . The result is 
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and 
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The contribution A/4 arises from the terms contained in the structure functions W4 and W5 
in (78). Using the analytic results for the master integrals yields 
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Comparing the result for F P (A P ) in (93) with the expression for F E (A E ) in (92) we observe 
that, as anticipated, the cut on hadronic P + contains a much larger fraction of all B — > X u l~u 
events. In fact, Fp(Ap) is directly given in terms of an integral over the shape function, 
without a weight function of order Aqcd/^6- If we neglect radiative corrections for a moment, 
we simply have Fp(Ap) = f^ p duj S(u), where the shape function is expected to peak at a 
position Co m A rs 0.5 GeV. For the "optimal cut" P + < M^/Mp ~ 0.66 GeV, which eliminates 
the charm background entirely, we thus expect that more than half of all events are contained 
in the event fraction Fp(Ap). The hadronic P + spectrum therefore offers a very promising 
new avenue for a high-precision measurement of \V u b\. 

In a realistic measurement, it is often necessary to implement a loose cut on the lepton 
energy in order to eliminate soft leptons, which can be difficult to measure in the detector. 
Let us discuss how our results would change in the presence of such a cut. Assume that we 
require Ei > Eq with E$ of order 1 GeV or so, corresponding to x < xq = 1 — 2Eo/rrib with 
x = 0(1). Our expression for the event fraction Fp(Ap) in (93) remains valid, except that 
the master functions now become functions of the parameter x . We find that T(a), f n (a) and 
A/ 4 (a) in (94) and (95) get replaced by 

T{a, x ) = 2x 2 (3 + 2x ) [h{-a, x ) -h(-a, 1)] - 12x (l + x ) [h(l - a, x ) ~h(l- a, 1)] 



+ 12x 



I 1 (2-a,x )-h(2-a,l) + 6h(2 - a, x ) - 4^(3 - a, x ) , 



f n (a,x ) = 
Af 4 (a,x ) = 



[expression for T(a, xq) with l\ — > I n ] 



4x1 



T(a,x ) 
h(-a,x ) - I 4 (-a, 1) 



(97) 



6x 2 



I 4(1 — a, xo) — h(l — a, 1) + 2/ 4 (3 — a, xo) 
T(a,x ) 
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For x — 1 these complicated results reduce to the simpler expressions given above. Numer- 
ically, the cut on the lepton energy has a minor effect provided that Xq is larger than about 
0.6, corresponding to a lower cutoff E less than about 1 GeV. 



7.3 Hadronic invariant mass spectrum 

A cut y/sH < Mo on the hadronic invariant mass in the final state constitutes the ideal 
separator between B — » X u l~v and B — > X c l~v events, since any final state containing a 
charm hadron has invariant mass above Mb- The left-hand plot in Figure 3 shows that such a 
cut fully contains the region with P + < M D /M B considered earlier. In addition, the hadronic 
invariant mass cut contains a triangle-shaped region of larger P + , which culminates in a cusp 
where P + = P_ = Mp- Near the cusp, both light-cone momentum components are of the 
same order, and hence this portion of phase space lies outside the shape-function region. In 
other words, the region near the cusp is a dangerous one (a "Bermuda triangle"), where the 
theoretical description based on the collinear expansion breaks down. A priori, then, it is not 
evident that we can compute the fractional rate Fm(so) = F(sh < so)/r to t in a controlled 
heavy-quark expansion. 

To see what happens, it is instructive to first ignore radiative corrections. At tree level, it 
is straightforward to obtain 

F M (s ) = JdwS(cu)+ JduS{u)(^j ( 2 ~^)' (98) 

A s 

where A s = s /M B . The calculation of this event fraction requires knowledge of the shape 
function over a wider range in Co than in the case of the event fraction with a cut on P + . It 
is therefore more difficult to extract the relevant hadronic information from the B — > X s 7 
photon spectrum. The first integral is the same as in (93) and corresponds to the region in 
phase space where P + < s /M B . Note that the ratio A s = s /M B plays the same role as the 
cutoff Ap on the P + spectrum. The second integral corresponds to the phase space above the 
dashed line in Figure 3. The region near the cusp corresponds to the upper integration region 
in the second integral. Note that, due to the rapid fall-off of the integrand, the "Bermuda 
triangle" only gives a power-suppressed contribution to the decay rate and so can be ignored at 
leading order in the heavy-quark expansion. We will address this point in more detail below. 

When radiative corrections are included, the result for the integrated hadronic invariant 
mass spectrum becomes rather complicated. It is convenient to split up the integration region 
in phase space into a box-shaped region with P + < s /M B and a triangular shaped region 
with Sq/M b < P + < y/sQ. Writing 

Fm(s ) = F h M °*(s ) + F^ lc (s ) , with F h ™(s ) = F P (A S ) , (99) 

we find that the box contribution is given by the expression for the rate fraction Fp(Ap) in 
(93) evaluated with A P = A s = s /M B . For the remaining contribution from the triangular 
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region, we obtain 



j-i triangle , 
M 1 
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G 1 (AJoj) + ^^G 2 (A s ,oo) 
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(100) 



where 

Gi(z) = T(a, z){l+ CFa °[ mb) H (a, z) + ^E^M. [ 7 _ ^ _ 3^ z) + ^ z)] 1 (1Q1) 



An ' An 

contains the same functions T, H and /„ as defined in (94) and (95), but with all master 
integrals replaced by I n (b, 1) — > I n (b, z). In addition, we need 
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(102) 



6 Ji (2 - a, z) - Ah (3 - a, 2) = 4T(a, z) , 



fc 2 (z) = 4 6/2(2 - a, z) - 4/2(3 - a, z) - 3 6/1 (2 - a, z) - 4/1 (3 - a, 2) 



(103) 



As mentioned above, the phase-space region near the cusp where Cj ~ y^so or P ~ ^/so 
gives a power-suppressed contribution to the decay rate. Using the explicit results for the 
functions Gi together with the asymptotic form of the shape function in (77), we find that the 
corresponding term is 

F M (s ) 3 c ™™> ^iM _±_ (^) 3 - a (l + _A a 



7T (3 - a) 2 Vv/so/ V 4 3-a/ 
where the dots represent higher-order power corrections. Whereas the fraction Fm(so) is of 
0(1) in power counting, the result (104) scales like ( Aqcd/^) ( 3_a )/ 2 for s ~ ^Aqcd- in 
the derivation of the formula for the event fraction we have neglected other power-suppressed 
terms with the same scaling. For consistency, we should therefore omit the term in (104), 
which can be done by replacing all occurrences of ^/so in upper integration limits in (100) and 
(102) with 00. Only in that way we ensure that our calculations provide the unique leading- 
power contribution in the heavy-quark limit. We will use this prescription in our numerical 
analysis in Section 9. 

Because of the presence of power corrections from a region in phase space where the 
collinear expansion breaks down, it is not clear to us how one would construct a systematic 
heavy-quark expansion for the fraction Pm(so) beyond the leading order. Clarification of this 
point deserves further study. 
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Figure 4: Phase-space constraints (left) and weight functions (right) for com- 
bined cuts on the hadronic and leptonic invariant mass: (sq, g 2 ) = (M|>, 0) 
(solid), (M 2 D , 6 GeV 2 ) (dashed), and ((1.7 GeV) 2 , 8 GeV 2 ) (dotted). 

7.4 Combined cuts on hadronic and leptonic invariant mass 

Bauer et al. have proposed to reduce the sensitivity to shape-function effects in the extraction 
of \Vyh\ by combining a cut on hadronic invariant mass with a cut q 2 > q 2 on the invariant 
mass squared of the lepton pair [46]. The first plot in Figure 4 shows that this eliminates 
a large portion of the events with large P_. It is straightforward to study the effects of 
such a combined cut in the approximation where radiative corrections are neglected. For the 
corresponding event fraction at tree level, we obtain 

F comh (s H <s ,q 2 >q 2 )=yl(2-y ) J du S(u) + JduS(u)(^J ( 2 ~ |f ) ' (105) 

A s /yo 

where y = 1 — q 2 /(rribM B ), and A s = s /M B as above. For a fixed hadronic-mass cut s , the 
effect of the cut on q 2 is to broaden the support of the first integral, while at the same time 
reducing its weight due to the prefactor. To illustrate this point, we show in the second plot 
in Figure 4 the weight functions under the integral with the shape function for three different 
choices of (so, q^)- The sensitivity to the precise form of the shape function is reduced because 
the weight functions become progressively more shallow as the value of q 2 . is raised. However, 
this reduction comes at the price of a significant reduction of the rate, raising questions about 
the validity of the assumption of quark-hadron duality. We will see in Section 9 that the 
relative uncertainty due to shape-function effects is not strongly reduced when imposing an 
additional cut on q 2 . We are therefore not convinced that it is worth paying this price. 

7.5 Comparison with the literature 

Before concluding this section, let us comment on results for inclusive -B-decay spectra in 
the shape-function region previously published in the literature. Detailed predictions for 
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B — > X u l~P decay distributions including shape-function effects were presented in [47] and 
[18]. In these papers, 0(a s ) corrections were included at the level of the underlying parton 
spectra. In the last reference fully differential distributions are presented, which can be used 
to calculate arbitrary spectra and implement experimental cuts. Dedicated studies of the 
hadronic invariant mass spectrum can also be found in [48, 49, 50]. Similar analyses for the 
photon energy spectrum in B — > X s ^y decays were presented in [33, 51, 52]. In all these 
works, shape-function effects are implemented by convolving parton-model distributions with 
a primordial structure function, typically by replacing the 6-quark mass in expressions for the 
parton spectra by a new variable ml = m b + u [3]. This procedure is correct at tree level. 
However, it is no longer fully consistent when radiative corrections are included, because 
part of the 0(a s ) corrections in the expressions for parton-level spectra are absorbed into 
the renormalization of the shape function. In particular, this changes the sign of the leading 
Sudakov logarithm, as can be seen by comparing the coefficients of the lnp 2 /p 2 terms in (18) 
and (26). This point has also been emphasized in [14]. An attempt to include radiative 
corrections to the spectra in a systematic way was made in [53], where in particular the 
evolution of the shape function has been addressed. Since our evolution equation does not agree 
with the one found by this author [31], our results for decay spectra are also in disagreement. 

In [14] , an expression has been presented for the double differential rate in the variables Ei 
and Eh, which (apart from a typo) is consistent with our findings. However, in this paper no 
distinction between a s (m&) and a s (/ij) has been made in the next-to-leading order corrections, 
which entails significant perturbative uncertainties. Also, the important question of a physical 
definition of the 6-quark mass has not been addressed. Whereas in our case all quantities are 
defined in a physical subtraction scheme and no reference to the 6-quark mass is left in the final 
expressions for decay distributions (except as arguments of running couplings), the results of 
[14] contain explicit reference to the 6-quark pole mass. 

8 Model-independent relations between spectra 

The decay spectra and event fractions discussed in the previous section are given in terms of 
weighted integrals of perturbative expressions with a non-perturbative shape function, which 
encodes hadronic physics related to the bound-state properties of the B meson. To use these 
formulae, one must take recourse to a model for the shape function, or (better) extract the 
shape function from a fit to experimental data. Alternatively, it is possible to derive model- 
independent relations between different decay distributions in which the shape function has 
been eliminated [1]. This makes use of the fact that at leading power in Aqcd /mb shape- 
function effects in inclusive decays to light hadronic final states are described by a single 
universal (i.e., process independent) function. The most promising strategy is to relate event 
fractions in semileptonic B — > X u l~v decays to a weighted integral over the B — > X s ^ photon 
spectrum, which at present provides the most direct access to the shape function. 

While it is straightforward to derive such relations at tree level, radiative corrections in- 
troduce non-trivial complications [11, 12, 13]. In the following, we illustrate with a concrete 
example how such shape-function independent relations can be derived systematically within 
our framework. Since our formalism has yet to be applied to the B — > X s ^ photon spectrum, 
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we will instead derive a relation between the charged-lepton energy spectrum and a weighted 
integral over the P + spectrum in B — > X u l~v decays. This relation serves as a prototype for 
all other shape-function independent relations between partially integrated decay rates. Note 
that the properties of the P + distribution are very similar to those of the B — > X s ^ photon 
spectrum in the sense that, at tree level, both spectra are directly given in terms of the shape 

function, e.g. (l/T tot ) dT/dP + = S(P+) + This connection, and its potential usefulness 

for an extraction of |V^|, has been noted earlier in [26]. 

Specifically, we wish to construct a perturbative weight function w(A,P + ) such that at 
leading power in AQ CD /m b 

r M B /2 d y rA -p 

/ d E l —= dP + w(A,P + ) — , A = M B -2E . (106) 

This relation is independent of the shape function and hence insensitive to hadronic physics. 
The construction of the weight function is straightforward order by order in perturbation 
theory. Using the results of the previous section, we find 

/A „ N 2(A-P+) 3(4 -a) f C F a s (m b ) , . , 

w(A, P + ) = 4t 7 w x S 1 + , hUa) 

v ' +; M B -P + 6 -a 2-a | 4tt v; 



C F a s (ni) 
An 



h 2 {a) In mb(A ~ P+) + h 3 (a) 



(107) 



where 
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Large logarithms of the form [a s ln(mt//Xj)] ri and a s [a s ln(mf,//Xj)] n are resummed exactly at 
next-to-leading order in renormalization-group improved perturbation theory. They enter the 
coefficient functions through the parameters a and c defined in (79). 

In the equations for decay rates and spectra presented in Section 7, the dependence of 
the perturbative coefficients on the intermediate matching scale jii cancels against the scale 
dependence of the renormalized shape function S(u, fii). In practice, it is difficult to trace this 
cancellation if a model for the shape function at a fixed scale is employed. On the contrary, 
in the present case the weight function is formally independent of the scale /ij, because there 
is nothing to cancel a potential \Xi dependence in (106). This fact can also be shown explicitly 
using the formulae given above. Expanding the resummed result for the weight function to 
first order in a s , we obtain the simple expression 
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Figure 5: Weight function w(A,P + ) entering the rate relation (106) for A = 
/Mb and three different choices of the intermediate scale, namely fii = 
1.5 GeV (solid), 2.0 GeV (dashed), and 1.0 GeV (dotted). The weight function 
is formally independent of /ij. 

in which the dependence on \ii has canceled. However, since this formula contains a large 
logarithm and the scale to be used in a s is undetermined, it should not be used for phe- 
nomenological applications. 

Figure 5 shows results for the resummed weight function in the case where A = M 2 D jM B ~ 
0.66 GeV. The three curves refer to different values of the scale The stability with respect 
to variations of the intermediate scale is very good except for the case of a very low scale 
(fa = 1 GeV), for which the convergence of perturbation theory is expected to be poor. 

The result (107) settles an old argument about the form of the weight function in relations 
such as (106). Leibovich et al. have presented a form for the weight function in a relation 
between the resummed B — > X s ^ and B — > X u l~u decay rates, in which the kinematic 
variables corresponding to P + and A above enter in a most complicated form [11]. Their result 
has the unattractive feature that the integral over the weight function contains a Landau-pole 
singularity, which must be avoided by introducing a cutoff on the P + integral on the right-hand 
side of (106). A simpler form of the weight function, which is equivalent to (109) and contains 
only a single logarithm of the kinematic variables, has been promoted in [12]. Our exact result 
(107) obtained after renormalization-group improvement retains a simple form with only a 
single logarithm. Note, in particular, that our result does not exhibit any unphysical Landau 
singularities. We believe the resolution of the discrepancy has to do with the choice of the 
intermediate matching scale fii. Since the weight function is formally independent of fii, we are 
free to choose any value fii ~ ly/m&AQCD- Taking the particular choice fif = mt,(A — P + ) would 
eliminate the logarithmic term in (107), however at the price of introducing a very complicated 
dependence on the kinematic variables A and P + via the dependence of the weight function 
on the quantity a = || ln[a s ( / Uj)/Q; s (mf ) )]. Also, with this choice the weight function would 
develop a Landau pole singularity at the point where the coupling ct s (^/ m &(A — P+)) gets 
strong, which happens for (A — P + ) ~ AQ CD /m b . All these unwanted features are avoided by 
using a fixed value for the intermediate scale, as we did in (107). 
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9 Numerical results 



We are now ready to study the implications of our analysis for phenomenology. We start by 
deriving the numerical values for the shape-function mass and kinetic energy including errors. 
We then present a model for the shape function which satisfies all theoretical constraints, 
and study its behavior under renormalizat ion-group evolution. Finally, we present numerical 
results for the various decay rates and spectra investigated in Section 7. 

Throughout this paper we use the two-loop running coupling constant in the MS scheme, 
normalized such that a s (M z ) = 0.119. We take m b = 4.65 GeV as the default value for the 
6-quark mass in the shape-function scheme (see below), and ^ = 1.5 GeV as the standard 
choice of the intermediate matching scale. This corresponds to setting /z 2 = mf,A had with a 
typical hadronic scale Ah a d ~ 0.5 GeV. Note that this choice eliminates the appearance of the 
6-quark mass from the arguments of the logarithmic terms in the expressions for the decay 
rates. The values of the strong coupling evaluated at these scales are a s (m b ) ~ 0.222 and 
a s {ni) ~ 0.375. The corresponding values of the perturbative parameters a and c defined in 
(79) are a ~ 0.335 and c ~ 0.614. Finally, the leading-order Sudakov factor in (78) takes the 
value e VH{ - mb ^ ~ 1.21. 

9.1 Shape-function mass and kinetic energy 

A value for the shape-function mass can be obtained by combining the relations (69) or 
(70) with existing predictions for the 6-quark mass in the relevant renormalization schemes. 
The potential-subtracted mass at the scale /if — 2 GeV has been determined from moments 
of the bb cross section and the mass of the T(IS') state [54]. Using the first relation in 
(70), the result of this paper implies mf (2 GeV, 2 GeV) = m^ s (2 GeV) = (4.59 ± 0.08) GeV. 
From a similar analysis the kinetic mass has been determined at the scale /j,f — 2 GeV to be 
mf n (l GeV) = (4.57 ± 0.06) GeV [41]. From the second relation in (70) it then follows that 
mf F (l GeV, 1 GeV) = (4.65 ± 0.06) GeV. Using relation (75) to compute the scale dependence 
of the shape-function mass, we obtain at the intermediate scale the values m| F = 
(4.61±0.08) GeV and mf F (fii,fii) = (4.65±0.06) GeV, respectively. Alternatively, we may use 
relation (69) in conjunction with an experimental determination of the 6-quark pole mass from 
moments of inclusive B — > X c l~v and B — > X s 7 decay spectra. Using the average value A po i c = 
(0.375±0.065) GeV obtained from [55, 56, 57, 58], we find mf^m) = (4.67±0.07) GeV. It 
is quite remarkable that these different determinations of the shape-function mass, which use 
rather different physics input, give highly consistent results. Combining them, we quote our 
default value for the shape-function mass at the intermediate scale Hi = 1.5 GeV as 

m b w (fx i , Hi) = (4.65 ± 0.07) GeV . (110) 

The corresponding A parameter is A(/ij,/ij) = (0.63 ± 0.07) GeV. 

A value of the kinetic-energy parameter in the shape-function scheme can be obtained 
from (71) or (72). Using the first relation combined with the experimental value — Ai = 
(0.25 ±0.06) GeV 2 [56, 57, 58] yields /i 2 (/ii, ^) = (0.271 ± 0.064) GeV 2 . Alternatively, we may 
use the result for the kinetic-energy parameter obtained in the kinetic scheme, [/z 2 (1 GeV)] km — 
(0.45 ± 0.10) GeV 2 [41], to get from (72) the value /i 2 (/ii,/i;) = (0.254 ± 0.107) GeV 2 . Again, 
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the two determinations are in very good agreement with each other. Combining them, we 
obtain 

Ulim^i) = (0.27 ±0.07) GeV 2 . (Ill) 
9.2 Model shape functions 

In our analysis of decay rates below, we will adopt a model for the shape function S(u,fj,i) 
at the intermediate scale. We stress, however, that ultimately the shape function could be 
extracted from a fit to the photon spectrum in B — > X s 7 decays. This would largely reduce 
the theoretical uncertainties in our predictions. 

For the purpose of illustration, we use a two-component ansatz for the shape function that 
is a generalization of the model employed in [18, 33]. The form we propose is 

*<<*,/.) = X ( X J -p(-„-J — ^ ^-a i 21 "— + 1 > < 112 > 

where A and b are model parameters, and A differs from the pole-scheme parameter A po i e by an 
amount of 0(a s (/i)). In the limit a s (fj) — > this function reduces to the one used in [18, 33]. 
The radiative tail ensures the correct leading asymptotic behavior of the shape function as 
displayed in (77). This in turn gives the correct power-like dependence of shape- function 
moments on the integration cutoff. In our model, this tail is glued onto a "primordial", 
exponential function such that the combined result is continuous. 

There are several non-trivial constraints on the parameters of the model. The normalization 
factor N is given by 
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(113) 



which is determined such that the integral over the shape function from Cj — to /j,f + 
A(fj,f,n) coincides with the first expression in (74) up to second-order power corrections and 
exponentially small terms of order e -At/ / A , which are negligible whenever /if is sufficiently large 
to trust our moment relations. By evaluating the first moment of the model shape function, 
we find that 

A = A pol , + ^W^ (11 4) 
to first order in a s . From (69) it then follows that 

A = A(^, fa) i) Cpa ^ l) ^ K(ja, + 51 MeV . (115) 

Finally, the model parameter b can be adjusted to reproduce a given value for the second 
moment of the shape function. 

Table 1 collects the parameters of the model shape functions at the intermediate scale 
Hi = 1.5 GeV corresponding to different values of A(//j,//j) and fi 2 (//j, /ij), as computed from 
the moment relations (74). These quantities are varied within their respective error ranges 
given in (110) and (111). Note that the A values in the table are in very good agreement with 
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Table 1: Parameters and moments of the model shape functions at the inter- 
mediate scale Hi. The running quantities mf F , A, and /i 2 are defined in the 
shape-function scheme and evaluated at /if — /i — Hi — 1.5 GeV. 



Model 


Lines 


mf [GeV] 


A [GeV] 


A [GeV 2 ] 


A [GeV] 


b 


SI 


Dotted 


4.72 


0.56 


0.20 


0.611 


2.84 


S2 








0.27 


0.617 


2.32 


S3 








0.34 


0.626 


1.92 


S4 


Solid 


4.65 


0.63 


0.20 


0.680 


3.57 


S5 








0.27 


0.685 


2.93 


S6 








0.34 


0.692 


2.45 


S7 


Dashed 


4.58 


0.70 


0.20 


0.751 


4.40 


S8 








0.27 


0.753 


3.61 


S9 








0.34 


0.759 


3.03 



those obtained from the relation (115). The left-hand (right-hand) plot in Figure 6 shows three 
models for the shape function obtained by varying the parameters A and /i 2 in a correlated 
(anti-correlated) way. In both cases, the solid, dashed, and dotted curves refer to different 
values of A, as indicated in the table. 

In Figure 7 we illustrate how the shape function behaves under renormalization-group evo- 
lution. The sharply peaked solid line shows our model function evaluated with A = 0.495 GeV 
and b = 3.0, which we use as an ansatz for the function S(u>, /io) at the low scale /iq — 1 GeV. 
For comparison, the dotted gray curve shows the default choice for the shape function adopted 
in [18, 33], which exhibits a very similar shape except for the missing radiative tail. The broad 
solid curve gives the shape function at the intermediate scale ^ = 1.5 GeV as obtained from 
the evolution equation (59). The barely visible dashed-dotted curve shows our default model 
for the function S(u,fii), which coincides with the solid line in the left-hand plot. The beau- 
tiful agreement of the two curves gives us confidence in the consistency of our models adopted 
for the shape function at the intermediate scale. 

9.3 Predictions for decay spectra and event fractions 

We are now ready to present our results for the decay spectra and partially integrated event 
fractions in B — > X u l~v decays. In order to illustrate the sensitivity to shape-function effects 
we use the nine shape functions SI through S9 in Table 1, thereby taking into account the 
full range of allowed values for the parameters A(/Xj,/Xj) and /^(/ii, /Xj). For each physical 
quantity we draw three bands corresponding to the three different values of A. The width of 
each band reflects the sensitivity to the variation of /i 2 . While one might in principle consider 
other functional forms for the shape function (which, however, must be consistent with the 
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Figure 6: Various models for the shape function at the intermediate scale 
\ii = 1.5 GeV, corresponding to different parameter settings in Table 1. Left: 
Functions SI, S5, S9 with "correlated" parameter variations. Right: Functions 
S3, S5, S7 with "anti-correlated" parameter variations. 
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Figure 7: Renormalizat ion-group evolution of a model shape function from a 
low scale fiQ (sharply peaked solid curve) to the intermediate scale fii (broad 
solid curve). See the text for an explanation of the other curves. 



asymptotic behavior as predicted by the operator product expansion), we believe that the 
variation of our results corresponding to the various models presents a realistic estimate of the 
shape-function sensitivity (see below). 

The following predictions for spectra and rate fractions refer to the leading term in the 
heavy-quark expansion. Using our formalism, power corrections can be computed systemati- 
cally to any order in AQ CD /m b , by extending the two-step matching QCD — ► SCET — > HQET 
to the appropriate order. These power corrections fall into several categories, including correc- 
tions from phase space, kinematic factors, and subleading shape functions [59] (see [44, 60, 61] 
for tree-level discussions of these effects on various B — > X u l~v spectra). Estimating the 
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Figure 8: Fraction of B — > X u l~v events with hadronic light-cone momentum 
P + < A P (left), and fraction of events with hadronic invariant mass sp < s 
(right). In each plot, the three bands correspond to the values A = 0.63 GeV 
(solid curves), 0.70 GeV (dashed curves), and 0.56 GeV (dotted curves). Their 
width reflects the sensitivity to the value of ^ varied in the range between 0.20 
and 0.34 GeV 2 . The arrow indicates the point at which the charm background 
starts. 

corrections from phase space alone, we find that they typically change the leading-order pre- 
dictions for partially integrated event fractions by about 10%. A more careful investigation of 
power corrections is left for future work. Finally, we note that our calculations would break 
down if the cuts on kinematic variables were taken to be too strict, because then the spectra 
would become dominated by hadronic resonance effects. Parametrically, this happens when 
the quantities Ap, A g , or Ap become of order Aq CD /Mp ~ 50MeV. 

In Figure 8 we show predictions for the fractions of all B — ► X u l~v events with hadronic 
light-cone momentum P + < Ap, and with hadronic invariant mass squared sp < s . Recall 
that, for Ap = A s = s^/Mp, the hadronic invariant mass fraction Fm differs from the fraction 
Fp by the contribution of the events in the triangular region above the dashed line in Figure 3. 
Comparing the two plots, we observe that this additional contribution is predicted to be very 
small. (Note that for large values of A s we even find a negative contribution to the rate from 
the triangle region for some choices of the shape function. This feature is unphysical and 
should be fixed by the inclusion of power corrections to our leading-order predictions.) The 
arrows on the horizontal axes indicate the points Ap s = Mp/Mp, beyond which final states 
containing charm hadrons are kinematically allowed. With this choice of the cut, both rate 
fractions capture about 80% of all events. While it is well known that a hadronic invariant mass 
cut y/sp < Mp, provides a very efficient discrimination against charm background [48, 49, 50], 
here we observe that the same is true for a cut on the P + variable. Cutting on P + offers 
the additional advantage of a "buffer zone" against charm background. Whereas the region 
in which charm final states are kinematically allowed borders the region with ^fsp < Md> if 
touches the phase-space region with P + < M^/Mp at only a single point (see Figure 3). 

Our results for the charged-lepton energy spectrum Si(Ei) = (l/T tot ) (dT/dEi), and for 
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Figure 9: Charged-lepton energy spectrum in the region near the kinematic 
endpoint (left), and fraction of events with charged-lepton energy E\ > E 
(right). The meaning of the bands and the arrow is the same as in Figure 8. 

the event fraction with a cut E\ > Eq, are displayed in Figure 9. The right-hand plot shows 
that with Ae = M^/Mb only about 10-15% of all events are retained, and the theoretical 
calculation is very sensitive to shape-function effects. Such a cut is therefore much less efficient 
than the cuts on hadronic invariant mass or P + . As a result, an extraction of \V u b\ from the 
charged-lepton endpoint region is theoretically disfavored. 

Let us now comment in more detail on the quantitative features of the results for the various 
event fractions. The most significant observation drawn from Figure 8 is that the shape- 
function sensitivity is rather small for values of Ap and A s near the charm threshold. This is 
to some extent a consequence of our improved knowledge of the shape-function parameters. 
For instance, whereas in previous analyses the heavy-quark mass was varied in the range 
mi, = (4.8 ±0.2) GeV [18], the mass defined in the shape-function scheme is known with much 
better accuracy, see (110). A similar statement applies to the width of the shape function and, 
more importantly, to its asymptotic behavior. Yet, the good convergence of the three bands 
in each plot seems puzzling at first sight, given that the model shape functions in Figure 6 
are still rather different at Co ~ 0.7 GeV. The event fraction Fp(Ap), in particular, is at tree 
level given by the area under the shape-function curves between and Ap. Remarkably, an 
interesting "focusing mechanism" arises beyond the tree approximation, which has its origin 
in a subtle interplay between the shape function and the jet function under the convolution 
integral in (93). The point is that the next-to-leading order corrections to the jet function 
(the terms proportional to a s (fii)) contain a double-logarithmic singularity at the endpoint of 
the integration domain (at Co = Ap), which comes with a positive coefficient. Consider now 
the integrals over the model shape functions from to Ap, assuming that Ap is well beyond 
the maximum of the curves. The function with the smallest (largest) area takes the largest 
(smallest) value at the endpoint. When the shape functions are weighted with the jet function, 
the logarithmic spike at Co = Ap gives a contribution to the integral that is proportional to 
S"(Ap) and so is largest (smallest) for the function with the smallest (largest) area. The net 
result is to balance the differences in the areas and make the integrals converge more quickly 
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Figure 10: Left: Four examples of shape functions with identical normalization 
and first two moments, corresponding to A(/Xj, fa) = 0.63 GeV and fa\(fa, fa) = 
0.27 GeV 2 , but different functional form. Right: Corresponding results for the 
event fraction Fp(Ap). 



toward a single value. The workings of this mechanism are nicely illustrated by comparing the 
results for the event fractions Fp and Fp in Figures 8 and 9. Whereas focusing takes place 
in the former case, the weight function under the integral for Fp in (92) vanishes at Cj — Ap, 
thereby suppressing the contribution from the logarithms in the jet function. This explains 
why no focusing is observed for the event fraction with a cut on charged-lepton energy. 

Another way of thinking about this mechanism is to notice that the broadening of the shape 
function under renormalization-group evolution from a low scale up to the intermediate scale 
(right-hand plot in Figure 6) is a perturbative effect, which should not lead to an increased 
shape-function sensitivity. Because the convolution of the shape function with the jet and 
hard functions is independent of the scale fa, the broadening of the shape function must be 
compensated by perturbative logarithms in the jet function. 

Note that it is crucial for this mechanism that the logarithmic terms in the jet function give 
a positive contribution near Cj = Ap. This focusing effect did not take place in earlier studies 
such as [18, 47, 48, 49, 50], where parton-model spectra were convoluted with a primordial 
shape function. As mentioned earlier, in the parton model the leading Sudakov logarithm 
comes with the opposite (negative) sign, hence causing an anti-focusing effect of the radiative 
corrections. This also explains why our prediction for the hadronic invariant mass fraction 
Fm exhibits a smaller shape-function sensitivity than what has been found in most previous 
analyses. 

In addition to their dependence on the moment parameters A and /i 2 , our predictions are 
also sensitive to the functional form adopted for the shape function. In order to study this 
sensitivity we have constructed four shape functions with identical zeroth, first, and second 
moments, but rather different functional forms. They are shown in the left-hand plot in 
Figure 10. Some of these choices are intentionally rather extreme, given what is known about 
the shape function from the measurement of the B — > X s 7 photon spectrum [55]. The point is 
to illustrate that even drastic variations of the functional form do not invalidate the predictions 
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Table 2: Comparison of different theoretical methods using inclusive 5-decay 
rates to extract the CKM matrix element \V u b\. The error on the efficiency 
represents the sensitivity to the shape function only All results refer to the 
leading term in the heavy-quark expansion. 



Method 


Cut 


Efficiency 


Hadronic invariant mass 


s H <M 2 D 
s H < (1.7 GeV) 2 
s H < (1.55 GeV) 2 


(81.413-2)% 

(78.2lt 2 9 )% 
(72.71°;*)% 


Hadronic P + 


P+< £| = 0.66 GeV 
P+ < 0.55 GeV 


(79.6l| 2 )% 
(69.0111)% 


Charged-lepton energy 


E i > = 2 - 31 GeV 
Ei > 2.2 GeV 


(12.5l|;|)% 
(22.2l| 2 )% 


Combined (sh, q 2 ) cuts 
[tree level only] 


s H <M 2 D , q 2 >0 
s H < M 2 , q 2 > 6 GeV 2 
s H < (1.7 GeV) 2 , q 2 > 8 GeV 2 


(74.6+1;!)% 
(45.7l^; 8 )% 
(33.4l}|)% 



presented earlier in this section. The right-hand plot in Figure 10 shows, as an example, the 
results for the event fraction Fp(Ap) obtained using the four functions. Comparing this plot 
with the corresponding one in Figure 8 shows that at present the variation of the parameters 
A and yU 2 covers the uncertainty in the functional form. We have checked that this is also true 
for the other decay distributions studied here. However, if in the future the errors on A and 
fi 2 were reduced significantly, one should also include the sensitivity to the functional form in 
the error estimate. 

A summary of our phenomenological results can be found in Table 2. For a variety of 
different experimental cuts we report the fractions of the contained B — > X u l~u events and 
indicate how these fractions vary under the variation of the shape-function models. In all 
cases the goal is to reject the charm background as efficiently as possible, while preserving 
a sufficiently large fraction of the signal events so as to obtain a reliable determination of 
the CKM matrix element \V ub \. We emphasize that our results refer to the leading term in 
the heavy-quark expansion, and that the rates are expected to be modified somewhat by 
power corrections. The uncertainties quoted on the contained event fractions reflect their 
sensitivity to shape-function effects only. No other theoretical uncertainties are included in 
these estimates. 

The first portion of the table contains results for the event fractions with a cut on hadronic 
invariant mass. The first line corresponds to the "ideal" cut < M D , which in theory elim- 
inates the charm background entirely. This cut retains about 80% of all B — > X u l~v events. 
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In practice, spill-over from final states containing charm hadrons in a realistic measurement 
requires to lower the cut on hadronic invariant mass to values slightly below M^. Two typical 
choices are also covered in the table. Lowering the cut reduces the contained event fraction 
by modest amounts, while the sensitivity to shape-function effects increases markedly. Yet, 
even with a cut at 1.55 GeV the uncertainty on |T4b| would still amount to only 4%. 

The second portion of the table shows two examples of cuts on the hadronic light-cone 
momentum P + . With the "ideal" cut P + < M^/Mb about 80% of all events are contained, 
which is as good as with the hadronic invariant mass cut. The sensitivity to shape-function 
effects is somewhat more pronounced but still at a very acceptable level. The corresponding 
uncertainty on | V u b\ is about 5%. Since the P + spectrum is directly related to the photon energy 
spectrum in B — > X s ^ decays, this uncertainty can be reduced further by using experimental 
information on the photon spectrum. Due to the fact that there is a "buffer zone" separating 
the signal events from the charm background, we expect less spill-over than in the case of a 
cut on hadronic invariant mass. It may therefore be experimentally feasible to implement the 
cutoff near the optimal value. However, even if the value of Ap must be lowered by a small 
amount, the efficiency remains high and the shape function sensitivity at an acceptable level 
(see the second entry in the table). 

The third portion of the table shows results for the case of a cut on the energy of the 
charged lepton. The efficiency is obviously much reduced in this case, even if the cut can 
be relaxed into the region where some charm background is present (second entry). The 
theoretical calculations are therefore more prone to uncertainties from other effects such as 
weak annihilation [62]. 

In the lower portion of the table we give results for some combined cuts on hadronic and 
leptonic invariant mass, which have been briefly discussed in Section 7.4. Contrary to the other 
cases, these numbers refer to the tree-level approximation and so should be taken with caution. 
For reference, we quote again the (tree-level) result for the pure hadronic invariant mass cut, 
which differs significantly from the corresponding result including radiative corrections. While 
the additional cut on leptonic q 2 reduces the shape-function sensitivity, it comes along with 
a strong reduction of the efficiency. For instance, the combined cut y/sji < 1.7 GeV and 
q 2 > 8 GeV 2 employed in a recent analysis of the Belle Collaboration [63] has an efficiency 
of about 33% (at tree level and leading order in Aqcd/^), which is much smaller than the 
efficiency of the pure hadronic invariant mass cut tJ~s~h < 1.7 GeV. However, the sensitivity to 
shape-function effects is only slightly better in the case of the combined cut. 

As a final remark, let us comment on the applicability of the theoretical framework de- 
veloped in this work. According to Figure 3 most of the B — > X u l~u events are located 
in the shape-function region of large P_ and small to moderate P + . Our approach is based 
on a systematic heavy-quark expansion valid in that region of phase space. It allows us to 
calculate inclusive decay rates integrated over domains AP_ ~ Mb and AP + <C Mb, where 
typically AP + ~ Aqcd- (In the examples above, AP + = Ae, Ap, or A s , respectively.) While 
the corresponding predictions for decay spectra and event fractions are sufficient to analyse 
experimental data over most of the phase space relevant to measurements of the CKM matrix 
element \V u b\, it would be of interest to extend the validity of the theoretical description out- 
side the shape-function region. In the case where AP_ ~ AP + ~ Mp are both large, the decay 
spectra can be computed using a local operator product expansion. An interesting question is 
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whether it will be possible to match these two approaches in some intermediate region of AP + 
values that are numerically (but not parametrically) large compared with Aqcd- If the two 
predictions were to agree in an overlap region, this could be used to construct a theoretical 
description of inclusive B — > X u l~v decay distributions that is valid over the entire phase 
space. While this is an exciting prospect, we note that performing a systematic operator 
product expansion in the overlap region is far from trivial. Because of the hierarchy of scales 
API <C AP + M B <C M|, again a two-step procedure is in order. An example of such an 
approach can be found in [64]. 

10 Conclusions 

We have calculated differential spectra and partially integrated event fractions for the inclusive 
semileptonic decays B — > X u l~u at next-to-leading order in renormalization-group improved 
perturbation theory, and at leading power in the heavy-quark expansion. The hadronic tensor 
entering the differential decay rates has been factorized into perturbatively calculable hard 
and jet functions, and J, which we give to one-loop order, and a universal shape function 
S containing non-perturbative physics below an intermediate scale fa ~ y^m^A^cD- This 
factorization has been obtained by matching QCD onto soft-collinear effective theory to inte- 
grate out hard fluctuations of order mj, and by matching the result further onto heavy-quark 
effective theory (HQET) to integrate out hard-collinear modes at the intermediate mass scale. 
Large logarithms have been resummed by solving the renormalization-group equation for the 
hard kernels, evolving the functions Hij from the hard scale fa ~ nib down to the intermediate 
scale fa. In order to make use of the resulting expression for the hadronic tensor, the shape 
function is needed at this intermediate scale fa. We have derived the anomalous dimension 
of the shape function at one-loop order and given an exact analytic solution of the result- 
ing renormalization-group equation. Our solution can be applied to any model for the shape 
function obtained at any scale. Alternatively, if the shape function is treated as a phenomeno- 
logical quantity, renormalizing it at the scale fa avoids doing perturbation theory at a low 
hadronic scale. 

Moments of the shape function are often identified with HQET parameters such as A and 
Ai. This identification was so far only understood at tree level. We have discovered the 
appearance of a radiative tail of the shape function that vanishes slower than This 
feature renders all shape-function moments (including the normalization integral) ultra-violet 
divergent. To define the moments consistently we have introduced a hard ultra-violet cutoff 
Auv on the integrals over uo. This is natural, since in any physical process the shape function 
will only be integrated over a finite interval. The dependence of the moments on this cutoff can 
be controlled using a local operator product expansion, which has enabled us to derive reliable 
perturbative relations between shape-function moments and HQET parameters. We have also 
obtained a formula for the normalization integral over the shape function as a function of 
Auv, and from it a model-independent prediction for the asymptotic behavior of the shape 
function for large values of |a>|. Surprisingly, this analysis reveals that the shape function has 
a negative tail and so is not positive definite, contrary to common expectation. 

Due to the fact that the 6-quark pole mass, and with it the HQET parameter A, suffers 
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from infra-red renormalon ambiguities, it is favorable to replace m\° 6 by a short-distance mass 
defined in a physical subtraction scheme. We have argued that the most natural definition 
for applications to inclusive spectra is to define a "shape-function mass" rrVjf (Auv , /•*) by 
enforcing that the first moment of the renormalized shape function vanish for any given value 
of the cutoff A uv . The dependence of the shape-function mass on the cutoff A uv and on the 
dimensional regularization scale \i is controlled by evolution equations, which can be trusted 
as long as both scales are much larger than Aqcd- Using a similar approach, we have defined 
a running kinetic-energy parameter /^(AuvjZ-O m the shape-function scheme, which can be 
used to replace the (ambiguous) HQET parameter Ai. We relate our new parameters to some 
previous, physical definitions of m b and //^. All of these new insights restrict model building 
of the shape function drastically. We give expressions for some model shape functions that 
are fully consistent with all constraints, albeit leaving enough freedom to accommodate future 
experimental constraints from the B — > photon spectrum. 

In the second part of the paper we have applied our results to make predictions for several 
interesting B — > X u l~u decay distributions and spectra in the shape-function region of large 
hadronic energy {Eh ~ Mb) and small hadronic invariant mass {sh ~ M#Aqcd)- We have 
presented results for event fractions with cuts on the charged-lepton energy, hadronic invariant 
mass, or hadronic light-cone momentum P + = Eh — \Ph\, where in each case the cuts are 
chosen so as to reject background from B — > X c l~v decays. These formulae are presented 
as convolution integrals of weight functions with the shape function renormalized at the in- 
termediate scale Our results are valid at next-to-leading order in renormalization-group 
improved perturbation theory and at leading power in the heavy-quark expansion. The dis- 
cussion of these distributions is most transparent in terms of the hadronic phase space for the 
two variables P± = Eh =F \Ph\, which we have studied in some detail (see Figure 3). 

An important outcome of our phenomenological analysis is the finding that a cut P + < 
M^/Mb — 0.66 GeV on the hadronic light-cone momentum eliminates the charm background 
while containing the vast majority of all B — > X u l~v events. At leading power in Aqco/^fc 
this cut has an efficiency of about 80%, which is almost as efficient as the cut y/sn < Mo 
on hadronic invariant mass. A cut on P + offers several advantages over a cut on hadronic 
invariant mass. First, it provides a "buffer zone" against charm background, which borders 
the signal region at only a single point in phase space. Secondly, the hadronic physics affecting 
the P + spectrum is directly related to the hadronic physics affecting the B — > X s 7 photon 
energy spectrum (the two spectra are identical at tree level). This implies that a simple, shape- 
function independent relation between the two distributions can be derived, which could be 
used to eliminate hadronic uncertainties (at least at leading power in Aqcd/^6)- Finally, as 
we have pointed out, the calculation of the hadronic invariant mass distribution suffers from 
the fact that it includes a region in phase space where the collinear expansion breaks down. 
While this region gives only a power-suppressed contribution to the rate, its presence might 
cause complications if the calculation is taken beyond the leading power. The P + spectrum, 
on the other hand, can be accurately calculated beyond the leading order using the methods 
developed here. It would be a most useful quantity to measure. 

Using our formalism, shape-function independent relations between different decay distri- 
butions can be derived in a systematic way, and they are free of unphysical Landau singulari- 
ties. As an example, we have derived a relation between the charged-lepton energy spectrum 
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and a weighted integral over the P + spectrum in B — > X u l~u decays. This formula serves as 
a prototype for other relations, including the more useful relation between the P + spectrum 
and the photon energy spectrum in B — > X s 7 decays, which can be used for a new, high- 
precision determination of \V u b\ once the technology developed in this paper has been applied 
to B — > X s 7 decays. 

Finally, let us reiterate that the results presented here are valid at leading power in 
Aqcd/"^?)- While it may be laborious to calculate the next term in the expansion, this can be 
systematically done in our formalism (with the possible exception of the hadronic invariant 
mass spectrum), and we believe it must be done in order to achieve a theoretical accuracy in 
the prediction of decay spectra at the 10% level. If this level of precision can be achieved, 
it would for the first time lead the way toward a precision measurement of the CKM matrix 
element \V u b\ with a theoretical error at the 5% level. This would be worth the effort. 

Note added: After this work was completed, we became aware of a paper by Grozin and 
Korchemsky [65], where the evolution kernel for the shape function in (45) was calculated at 
two-loop order. Our one-loop results in (46) agree with their findings. 

Acknowledgments: This research was supported by the National Science Foundation under 
Grant PHY-0098631. 
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